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1» Outline  of  Problem. 


One  of  the  striking  results  of  modern  astrophysics  is  the 
ubiquity  of  magnetic  fields  in  the  universe.  At  present  we  know 
a large  number  of  stars  with  magnetic  fields  of  the  order  of  sev- 
eral thousand  gauss  (Babcock  and  Cowling,  1955).  A rather  largo 
fraction  of  them,  if  not  all,  have  variable  fields,  so  much  so 
that  one  is  tempted  to  consider  variability  as  an  intrinsic  feature 
of  magnetic  stars.  Moreover,  there  is  evidence  that  the  clouds 
of  rarefied  gas  which  are  found  in  galactic  space  carry  magnetic 
fields.  The  most  successful  theories  ox  the  origin  of  cosmic 
rays  seem  to  be  those  that  assume  the  acceleration  of  cosmic- 
ray  particles  to  be  caused  by  the  mean  action  of  these  fields 
(Fermi,  1949 j Morrison,  Clbert,  and  Hnssi,  1954). 

Thu  theory  of  hydromagnetism,  or  magnetohydrodynamics, 
has  only  fairly  recently  been  developed  in  an  effort  to  explain 
these  cosmic  magnetic  fields,  and  it  bids  fair  to  provide  such  an 
explanation.  As  we  shall  see  (Sec.  4)  small  stray  magnetic  fields 
can  be  amplified  by  the  action  of  suitable  fluid  motions.  Now 
fluids  of  large  dimensions  are  as  a rule  highly  turbulent . One 
might  therefore  expect  the  velocity  distribution  of  those  fluids 
to  follow  ecmo  statistical  pattern,  so  that  the  amplification  of 
stray  magnetic  fields  may  be  determined  from  statistical  prin- 
ciples. The  final  result  of  such  randomly  distributed 

The  worua  hydromagriet ism  and  magnetohydrodynamics, 
and  the  corresponding  adjectives  have  been  used  rather  indis- 
criminately in  the  literature.  We  are  using  here  the  former 
term  for  aesthetic  reasons  and  as  being  more  economical  of  space. 
This  usage  has  the  approval  of  the  distinguished  Secretary  of 
the  American  Physical  Society. 


ramplif icatory  processes  would  then  be  a statistical  equilibrium 
between  the  turbulent  motion  and  the  more  or  less  irregular  mag- 
netic fields  generated.  This  is  the  approach  taken  by  Batchelor 
(1950)  who  succeeded  in  estimating  the  magnetic  spectrum  that 
would  be  in  equilibrium  with  the  velocity  spectrum  of  turbulence. 
Such  a picture  gives  us  a first  insight  into  the  mechanism  of 
the  generation  of  magnetic  fields  under  the  conditions  studied 
by  the  astrophysicist.  It  indicates  in  a provisional  fashion  at 
least,  that  the  hydromagnetic  theory  can  account  for  the  presence 
of  magnetic  fields  in  the  universe  without  ad  hots  assumptions} 
the  generation  of  magnetic  fields  of  the  order  of  magnitude  ob- 
served follows  without  difficulty  from  the  application  of 
Maxwell's  electromagnetic  field  equations  to  moving,  electrically 
conducting  fluids  of  large  dimensions  (Sec.  2). 

From  this  viewpoint  the  older  experiences  regarding 
magnetic  fields  of  tre  earth  and  the  sun  appear  in  a new  light. 
They  are  taken  out  of  their  conceptual  isolation  and  appear  ss 
special  cases,  relatively  more  accessible  to  our  observation,  of 
a universal  phenomenon.  The  pertinent  facts  concerning  the 
earth's  interior  which  form  the  physical  background  for  the 
earth's  magnetic  field  have  been  reviewed  in  detail  some  years 
ago  (Elsaaaer,  1950).  We  shall  confine  ourselves  here  to  a few 
words.  The  earth  has  been  shown  from  seismological  observations 
to  have  a liquid  core  set  off  from  its  outer,  solid  part  by  a 
3harp  surface  of  discontinuity.  The  radius  of  the  core  is  about 
5500  km,  corresponding  to  55#  of  the  earth's  radius.  Geochemical 
evidence  indicates  that  tho  material  constituting  the  core  is 
primarily  molten  iron  with  perhaps  an  admixture  of  nickel  and 
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possibly  some  minor  constituents  In  solution.  In  order 
hydromagn e tic  theory,  to  the  earth’s  core  it  suffices  to  assume 
that  the  material  is  a good  electrical  conductor*  its  conductivity 
being  comparable  at  least  in  order  of  magnitude  to  that  of  ordinary 
metals.  Beyond  this  it  is  merely  necessary  to  assume  that  the 
core  is  fluid  and  that  internal  motions  occur.  The  geomagnetic 
secular  variation  puts  the  fluid  character  of  the  core  in  evi- 
dence. The  secular  variation  may  be  analyzed  into  a spectrum 
whose  prime  components  have  periods  of  the  order  of  a few  hundred 
years.  There  is  no  known  way  of  accounting  for  periods  of  this 
order  on  the  basis  of  mechanical,  thermal  or  other  processes 
occuring  in  the  solid  outer  parts  of  the  earthy  whereas  the  theory 
which  assumes  this  secular  variation  to  be  associated  with  fluid 
motions  in  the  earth's  core  is  able  to  explain  them  quite  satis- 
factorily (Elsasser,  1950).  In  the  present  review  we  shall  not  yet 
discuss  these  phenomena  in  detail  (although  they  are  well  suited 
to  compare  the  theory  in  a quantitative  fashion  with  direct 
observations).  Instead,  we  shall  focus  our  attention  on  the  basic 
hydromagnetic  processes  by  which  the  earth's  dipole  field  is 
generated  and  maintained.  These  processes  are  net  necessarily 
the  ones  which  are  most  directly  revealed  by  the  secular  variation? 
they  take  place  In  the  deeper  parts  of  the  earth's  core,  whereas 
the  geomagnetic  secular  variation  may  be  shown  to  inform  us  only 
about  the  conditions  in  a very  shallow  layer  of  fluid  adjacent 
to  the  surface  of  fcha  core. 

The  driving  mechanism  by  which  the  fluid  motion,  and  hence 
Indirectly  the  magnetic  field.  Is  maintained  is  generally  assumed 
to  be  thermal  convection  (Bullard,  1949)  although  convective 
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motion  induced  by  other  means,  in  particular  progressive  sedi- 
mentation (Urey,  1952)  need  not  be  rejected.  In  any  event,  the 
power  supplied  from  thermal  sources  ( radioactive  heat  in  the  core, 
plausible  radial  temperature  gradients)  is,  even  under  very  con- 
servative assumptions , more  than  sufficient  to  maintain 
throughout  the  lifetime  of  the  earth  fluid  motions  of  the  mag- 
nitude inferred  from  the  observed  secular  variation  (0.1-1 
xnm/sec).  The  details  of  the  primary  driving  mechanism  need  not 
concern  us  in  this  review}  it  will  appear,  however  (Sec.  5)  that 
the  fluid  motions  must  be  essentially  three-dimensional;  a pattern 
restricted,  e.g.,  to  spherical  sheets  or  to  meridional  planes  is 
not  adequate  to  produce  dynamo  action.  Hecent  seismological  re- 
search (Bullen,  IS 54)  makes  it  likely  that  the  central  part  of 
the  core  is  again  solid,  but  t;he  volume  of  this  inner,  solid 
sphere  is  only  a very  small  fraction  or  the  entire  volume  of  the 
core,  and  for  the  purposes  of  the  analysis  given  below  it  will 
be  sufficient  to  assume  the  earth's  core  as  a homogeneous  fluid 
sphere.  Compressibility  effects  are  not  likely  to  be  important 
and  ao  the  fluid  dare  m ay  be  considered  as  incompressible)  moreover 
the  electrical  conductivity  will  be  assumed  constant. 

There  is  another  class  of  extensively  studied  phenomena 
which  can  be  attributed  to  hydromagnetic  effects,  namely,  the  mag- 
netic fields  observed  on  the  sun,  particularly  , in  sunspots. 

Every  sunspot  has  a magnetic  field  associated  with  itj  the  larger 
the  spots,  the  larger  as  a rule  the  fields.  The  field  strength 
in  the  larger  spots  goes  up  to  a saturation  valuo  of  about 
3000  gauss  with  a margin  of  fluctuation  of  nearly  + 1000  gauss.-  Sun- 
spots appear  very  frequently  in  pa~rs,  the  line  connecting  the 
two  spots  running  in  an  east-westerly  direction,  along  a circle 
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of  latitude.  The  two  members  of  a sunspot  pair  always  have  op- 
posite magnetic  polarity.  The  " leader”  spot  (the  one  appearing  a- 
head  in  the  sense  of  the  solar  rotation)  always  has  ono  definite  pol- 
arity during  one  and  the  same  sunspot  cycle,  the 'Tollower”  havipg 
the  opposite  polarity.  The  le  ader  is  in  the  average  larger  and  has  a 
larger  magnetic  field  than  the  follower;  single  spots  as  a rule  have 
the  same  magnetic  polarity  as  the  leader  spots.  Sunspots  appear  and 
grow  to  their  full  size  in  the  course  of  a few  daysj  they  then  gradu- 
ally decay  during  the  course  of  a few  weeks,  some  of  them  persisting 
over  a few  months.  The  solar  latitude  at  which  the  spots  are  seen  is 
a function  of  the  11  1/2 -year  sunspot  cycle.  In  the  beginning  qf 
the  cycle  the  spots  appear  at  a latitude  of  about  30°.  As  the  sun- 
spot cycle  progresses  the  spots  appear  at  lower  and  lower  latitudes 
until  towards  the  end  cf  the  cycle  they  are  found  very  near  the  solar 
equator,  at  latitudes  of  5-10°.  At  the  same  time  new  spots  begin  to 
appear  around  the  latitude  of  30°,  but  these  spots  have  opposite  mag- 
netic polarities.  Curing  the  entire  subsequent  cycle  the  polarity  of 
th8  spots  Is  the  opposite  of  that  found  in  the  previous  cycle;  for  the 
next  11  l/2-ye«r  cycle  the  polarity  reverses  again,  and  so  on.  Clearly 
the  complete  sunspot  cycle  must  contain  the  reversal  of  the  field* 
and  extends  over  23  years.  There  are  many  other  observations  tfyat. 
indicate  the  presence  of  magnetic  fields  in  the  sun  other  than 
sunspot  fields  and  the  general  character  of  these  fields  appear^ 
to  vary  with  the  sunspot  cycle,  but  these  fields  are  very  much 
smrller  in  magnitude  than  che  sunspot  fields.  We  may  refer  here 
to  a very  comprehensive  recent  work  on  the  sun  (Kuiper,  ed,  1953)* 

Let  us  now  return  to  our  initial  remark  that  magnetic 
fields  in  cosmic  fluids  may  be  produced  by  amplification  from 


small  initial  fields*  We  referred  to  turbulence  as  being  able  to 
generate  such  fields  in  a statistical  manner*  Now  while  the  earth’s 
field  is  highly  irregular  in  the  detail by  far  the  largest  part 
of  the  field  has  the  form  of  a dipole  roughly  parallel  to  the 
earth’s  axi3.  The  prominent  feature  of  solar  magnetism  is  the 
23-year  cycle  which,  while  subject  to  fluctuations  ir  strength 
and  also  to  certain  fluctuations  in  length,  has  been  observed  to 
occur  with  consistent  regularity  over  the  last  200  years,  and  there 
5s  no  reason  to  doubt  that  it  is  a relatively  stable  feature  of 
solar  activity.  Similarly,  when  stars  are  observed  to  have  over- 
all fields  of  several  thousand  g&U33,  a systematic  cause  mu3t  be 
operative.  It  is  clear  that  if  the  observed  phenomena  are  to  be 
explained  by  hydromagnetic  amplification,  some  regularity  of 
the  pattern  of  fluid  motion  must  underly  them.  Fluid  motions  which 
produce  relatively  stationary  or  periodic  magnetic  fields  will 
be  designated  as  hydromagnetic  dynamos.  Since  fluid  motions  in 
large  dimensions  are  turbulent, 'J7  at  least  more  cr  less  irregular, 
we  shall  not  be  concerned  with  rigorous  solutions  of  the  hydro- 
magnetic  equations,  but  with  typical  solutions  which  demonstrate 
the  stability  of  the  fields  in  the  mean.  This  point  of  view  is 
in  full  agreement  with  the  observations  which  show  that  none  of 
the  observed  parameters  of  the  field  i3  rigorously  constant; 
all  of  them  are  subject  to  certain  rather  appreciable  margins 
of  fluctuations.  As  an  example  of  this  we  may  mention  that  the 
earth’s  dipole  moment  ha3  decreased. by  about  since  1850 

(Elsasser,  1950)  although  other  observational  data  leave  no  doubt 
that  the  earth  must  have  possessed  a magnetic  dipole  moment  for  a 
very  long  time  indeed,  so  that  the  present  decline  is  in  all 


likelihood  only  a temporary  fluctuation.  A disregard  of  these 
pronouncedly  irregular  fluctuations  of  the  field  has  led  to  the 
so-called  "fundamental"  theories  which  try  to  relate  geomagnetism 
and  other  cosmic  magnetic  phenomena  to  properties  of  matter  in 
the  large  not  contained  in  the  conventional  equations  ox  classical 
physics.  The  hydromagnetic  theory  shows  that  classical  physics 
can  account  for  these  phenomena,  but  that  the  particular  para- 
meters or  combination  of  parameters  have,  for  good  and  sufficient 
reason,  escaped  observation  in  the  laboratory  (Sec.  2). 

The  dynamo  theory  requires  that  the  fluid  motions  exhibit 
certain  regularities  in  order  that  magnetic  fields  may  be  main- 
tained in  the  average.  In  other  lADrds,  there  must  be  some  order- 
ing principle  that  controls  the  fluid  motions,  and  we  must  identify 
this  principle.  One's  first  idea  would  seem  to  be  a search  for 
arguments  of  symmetry : to  find  some  symmetry  requirement  that 

restricts  the  generality  of  the  fluid  motions  and  impresses  upon 
them  a relatively  simple  pattern.  But  such  a search  proves  to  be 
in  vain.  There  is  every  indication  that  an  appreciable  degree 
of  symmetry  of  the  fluid  motions  will  suppress  or  cancel  the 
effects  of  hydromagnetic  amplification  (Sec.  5)j  thus  we  are 
led  to  look  for  patterns  of  the  fluid  motion  of  a low  rather  than 
a high  degree  of  symmetry.  To  make  a long  story  short,  the  order- 
ing principle  which  engenders  the  most  conspicuous  hydromagnetic 
effects  may  be  identified  as  the  Coriolis  force  acting  upon  the 
fluid  motiona  in  a rotating  system.  Such  is  the  working  hypothesis 
of  this  paper.  The  Coriolis  deflection  affects  the  fluid  motions 
in  such  a fashion  that  the  resulting  pattern  does  not  in  general 
admit  of  any  symmetry  operations.  Stationary  or  periodic 


hydromagnetic  amplification  is  thus  related  to  the  rotation  of 
the  fluid  mass  in  which  it  occurs.  There  is  some  observational 
presumption  in  favor  of  this  idea,  since  the  star.s  with  strong 
magnetic  fields  seam  to  rotate  rather  rapidly.  The  evidence  is, 
however,  not  entirely  conclusive  and  the  assumption  must  be  jus- 
tified by  working  out  its  dynamical  consequences j we  might  remark 
that  we  are  not  informed  of  any  other  dynamical  principle  which 
could  be  adduced  to  explain  magnetic  fields  of  the  type  observed. 
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2 * Field  Equations,,  Dimensions, 


We  shall  use  the  rationalized  raks  system  throughout*  The 


electromagnetic  fields  will  be  assumed  to  obey  Maxwell's  equations, 
hence  the  material  medium  may  at  any  point  be  described  by  the 
three  constants  <p,  ix,  £.  . We  whall  assume  |x  and  £ as  constant 
throughout  space  and  shall  for  simplicity  assume  a-  constant  for  a 
given  fluid,  although  it  would  not  bs  difficult  to  generalize  the 
theorv  to  fluids  with  variable  cr  • We  then  have 


V x £ - -aB/at  , V • e = nA 

B = 0 , 


(2.1) 
' (2.2) 


where  ^ are  charge  and  current  density.  Next  we  write  down 

the  most  general  expression  for  the  current  density  admissible  in 
Maxwell < a theory 

=<ris  + <rv  x b + e.3E/at  + ^v  (2.3) 

fV* 1 / » • ^ v' 

where  v is  the  material  velocity  of  the  fluid  and  where  the  terms 
on  the  right  represent,  respectively,  the  conduction  current,  the 
induction  current,  the  displacement  current,  and  the  convection 
current.  Formula  (2.-3)  differs  from  the  conventional  expression 
for  the  total  current  by  the  second  and  fourth  terms  on  the  right 
which  contain  the  fluid  velocity,  v.  For  the  detailed  derivation 
of  these  terms  any  extensive  text  on  electrodynamics  may  be 
consulted. 

We  shall  now  proceed  to  show  that  the  last  two  terms  of 
(2*3)  are  negligibly  small  under  the  conditions  met  with  in  cosmic 
fluids.  We  shall  use  braces,  ^ , to  designate  the  order  of 

magnitude  of  a given  physical  quantity;  let  in  particular  { ?»•  ) 
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stand  for  a typical  length  and  ^co  j for  a typical  reciprocal  lime. 

We  first  note  that 

{v/c)  •=  r p j « i (2.4) 

• si^e  the  velocities  of  cosmic  fluids  rarely  exceed  a few  km/sec. 

We  now  compare  the  displacement  current*  the'  third  term  on  the 
right-hand  side  of  (2.5),  to  the  conduction  current,  the  first 
tern.  The  ratio  is 

(c ofi/trj  = (Y  J (2.5) 

*7 

'where  / is  exceedingly  small.  To  shew  this,  let  <T  = 10',  the 

18 

conductivity  of  ordinary  iron}  then  for  Y = 1 we  find  to  #s-10  . 

This  shows  that  for  the  frequencies  of  all  macroscopic  motions 
Y is  utterly  negligible.  (The  quantity  (2.5)  is  familiar  to  the 
student  of  metal  optics  where  it  is  used  in  the  same  way  as  here, 

« namely,  to  measure  the  ratio  of  displacement  to  conduction  current. ) 

It  is  readily  shown  that  the  ratio  of  convection  current  to  con- 
auction  current  is  also  given 'by  (2.5).  Prom  (2*1)  re  have  indeed 
1^7)  = ^£.EA}  j hence  this  ratio  is 

|^v/<rEj  = = {Y  j (2.6) 

i 

if  we  identify  w = v/X  as  a typical  frequency  of  the  material 
motion  of  the  fluid. 

We  next  find  for  the  ratio  of  the  electrical  to  the 
magnetic  field  energy,  using  (2.1) 

f£E2/V'1B2)  = ^E2/c2B2  j = [x2»2/c2J  = (p2J  (2.7) 

if  we  again  identify  X&>  with  the  velocity  of  the  fluid,  as  is 
proper  in  an  entirely  Maxwellian  scheme.  Hence  the  electrostatic 
field  energy  is  small  and  it  follows  from  familiar  arguments  that 
in  our  approximation  all  electromagnetic  processes  are  aperiodic. 


These  estimates  call  for  some  further  comments.  We  see 
from  the  numerical  estimate  of  Y that  even  for  a moderately  ion- 
ized gas  the  electromagnetic  phenomena  are  aperiodic  for  frequencies 
in  the  radio  spectrum*  Thus  plasma  oscillations  at  these  fre- 
quencies require  that  the  description  in  terms  of  the  macroscopic 
equations  of  Maxwell's  theory  be  invalid}  this  is  because  (2.7) 
becomes  invalid  for  the  velocity  of  the  electronic  component  of 
the  plasma.  In  the  present  article  we  are  dealing,  however,  cnly 
with  the  macroscopic,  average  motions  of  the  conducting  material* 
the  characteristic  frequencies  are  then  lower  by  many  powers  of 
ten  than  the  frequencies  of  the  radio  spectrum*  hence  (2.5), 

(2.6)  and  (2.7)  are  certainly  small  and  our  approximation  may  safely 
be  applied.  Now  (2.2)  and  (2.3)  give 

(2*3) 

Continuing  our  dimensional  analysis  (Elsasser,  1S54)  we 
compare  the  order  of  magnitude  of  the  three  terms  in  (2.8).  We 
first  notice  that  by  virtue  of  (2.7)  the  two  terms  on  the  right 
hand  side  are  of  comparable  order  of  magnitude,  -he  ratio  of  any 
one  of  these  terms  to  the  net  current  on  the  left  is 

r^TXv]  - (Rm\  (2.9) 

where  tne  non-dimensional  quantity  Rm  will  be  designated  as  the 

rwocrnftM  c Reynolds  number.  If  we  substitute  numerical  values  for 

the  quantities  on  the  left  of  (2.9)  we  find  that  R is  numerically 

m 

large  for  cosmic  fluids.  It  is  of  the  order  of  magnitude  of  sev- 
eral hundred  to  perhaps  a thousand  for*  the  earth’s  metallic  core, 

5 

depending  on  the  detailed  assumptions  made,  and  is  of  order  10 
or  more  for  most  astrophys ical  conditions.  This  constitutes  the 


Ij 
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essentlal  difference  of  cosmic  hydromagnetism  from  laboratory  con- 
ditions where,  as  one  readily  verifies,  is  numerically  small. 
Hence 


E — v a B 


(2 •10) 


Ui  r-t  o V. 


xum  •»  xx 7 


i Vk-rr  +*  V» 


•1‘  *-*  *»  ux 


ca  magnetic  field  is  maintained  is 


and  t ilS  luG  c han 

quite  at  variance  from  the  conventional  situation  where  the  net 
current  is  the  source  of  B. 

Operating  with  the  curl  on  (2.8)  and  using  the  field 
equations  we  eliminate  and  find 


SB /at  = X7x  (v  x B)  + v„V2B 

Ml  — IH  — 


(2.11) 


where  we  have  written 


v = (tier)"1  (2.12) 

m 

The  quantity  v will  be  designated  as  the  magnetic  viscosity. 

in  — — — ■ 

We  see  from  (2.9)  that  Rm  differs  from  the  conventional  hydro- 
dynamic  Reynolds  number,  R,  only  in  that  replaces  the  kina- 
matical  viscosity,  v,  of  the  fluid* 

The  physical  implications  of  (2.11)  are  best  brought  out 
in  terms  of  an  integral  equation.  To  obtain  it  we  integrate  (2.1) 
along  a contour  C and  use  (2.8)?  then  by  Stokes ? theorem 

(e/St)/Bdcr  = -J%  ♦ dC_ 


= B)  • dC_  - xj3)  • dC 

Now  if  the  first  integrand  on  the  right-hand  side  is  written 
B * (dC.xjr),  the  integral  can  be  given  a simple  geometrical 
meaning;  it  becomes  -/"B^dd”  where  the  integration  extends  over 
the  strip  that  the  contour  C subtends  in  its  motion  during  the 
time  dfc.  Since  y'fe^a  = 0 for  any  closed  surface,  we  can  write 


« 
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this 


( d/dt = -Vj/(^7  x B)  • d£  (2.13) 


On  the  left  there  appears  now  the  substantial  derivative,  refer- 
ring to  motion  with  the  fluid  particles.  Clearly,  if  (2.13)  is 
applied  to  any  arbitrary  contour,  it  is  equivalent  to  tho  dif- 
ferential equation  (2.11). 

Next  take  the  ratio  of  the  left-hand  side  of  (2.13)  to 
the  right-hand  side.  This  ratio  is  readily  sso  n from  (2.12)  and 
(2.9)  to  be  just  (r_|.  Hence  under  geophysical  and  astrophysical 
conditions  we  have  very  approximately 


(d/dt)^nd<r  = 0 (2.14) 

which  is  usually  enunciated  by  stating  that  the  magnetic  lines  of 
force  are  carried  along  bodily  with  the  fluid?  they  are  "frozen” 
as  it  were,  in  the  conducting  fluid.  On  applying  a well-known 
vector  identity  to  the  first  term  on  the  right  of  (2.11)  we 
Obtain  the  differential  equation  in  the  form 

dfi/dt  = (B  • X7  )v  - B(V  • v)  ♦ v V?'B 

which  may  be  further  simplified  on  introducing  from  the  equation 
of  continuity 

V7  • v_  = / d(f-1)/dt 

with  the  result 

d(f  _1B  )/dt  = Cf  • ^7)v  + r\  V2b  (2.15) 

^ JJJ  — +• 

an  equation  that  exhibits  more  clearly  the  role  of  compressibility 
(True 3dc 11 , 1950). 

This  equation  is  remarkable  in  that  it  shows  a complete 
formal  analogy  to  the  Helmholtz  theorem  of  the  conservation  of 
vortlcity.  Indeed,  if  we  replace  ES^  by  the  vorticity  vector  and 
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Vm  by  v,  (2.1b)  becomes  just  the  general  vortio ity-conservation 
theorem.  The  physical  implica lions  of  this  result  for  the  tubes 
of  vorticity,  which  now  may  be  transferred  at  once  to  the  tubes  of 
magnetic  flux,  are  well  enough  known  and  may  be  found  in  almost 
any  text  on  hydrodynamics.  There  is  one  point,  however,  where  a 
great  deal  of  misunderstanding  appears  to  exist  in  the  literature: 
One  has  become  habituated  to  saying  that  not  only  must  flux  be 
conserved  during  th9  motion  (for  vanishing  v or  vm)  but  that  also 
the'  lines  of  vorticitv  or  of  the  magnetic  field  must  be  closed 
in  the  absence  of  sources.  From  the  existing  discussions  one  is 
often  led  to  the  implicit  belief  that  the  condition  V • Br%  = 0 
requires  that  the  magnetic  lines  of  force  must  be  closed  curves. 
This  statement  is  certainly  incorrect.  The  subject  has  recently 
been  studied  in  some  detail  by  McDonald  (1954).  He  shows  that 
there  are  two  conditions  under  which  the  field  lines  are  not 
closed.  In  the  first  place  they  con  terminate  in  singularities 
that  is  points,  lines,  or  surfaces  where  J3  = 0.  Examples  of  or- 
dinary current  conf igurat ions  where  such  singularities  appear  can 
be  constructed  in  abundance}  similarly  such  singularities  may  be 
present  in  problems. of  hydromagnetism.  A second  class  of  non- 
closed  line3  is  that ' of  lines  which  are  "ergodic",  that  is  cover 
a region  everywhere  densely.  As  an  example  for  the  latter,  con- 
sider a current  system  consisting  of  a straight  current-carrying 
wire  and  a second  wire  forming  a circular  loop  which  lies  in  a 
plane  normal  to  the  straight  wir6  with  its  center  on  the  latter. 
The  lines  of  force  in  the  neighborhood  of  the  loop  are  clearly 
spirals}  a line  going  through  a given  point  will  be  closed  only 
when  the  ratio  of  the  currents  flowing  in  the  two  wires  has 
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certain  rational  values,  otherwise  the  line  will  fill  a torus- 
shaped surface  everywhere  densely.  It  cannot  be  our  aim  to  con- 
sider in  detail  these  somewhat  involved  properties  of  field  lines 
which  pertain  to  analytical  vector  geometry  and  are  of  more  in- 
terest perhaps  to  the  topologist  than  to  the  physicist*  one  may 
readily  p6~coive,  even  without  delving  into  more  rigorous  mathe- 
matics, that  the  closed  lines  of  force  form,  se t-theoretically 
speaking,  only  a subset  of  measure  zero  of  the  set  of  all  possible 
field  lines. 

One  thing  stands  out  clearly:  ^Intuitive”  arguments 

regarding  the  existence  of  closed  lines  of  force  and  the  impos- 
sibility of  generating  new  closed  lines  by  deformation  of  the 
fluid  are  of  no  value.  They  cannot  be  used  as  arguments  in  a 
discussion  of  hydromagnetic  amplification  unless  they  can  be 
converted  into  formulas  based  on  the  vector-field  equations  of  the 
theory.  Much  unjustified  scepticism  against  the  reality  of  hydro- 
magnetic  processes  has  arisen  from  such  intuitive  reasoning.  It 
can  be  avoided  only  on  abandoning  the  line -of -force  concept  of 
elementary  textbooks  in  favor  of  the  theory,  at  the  same  time  more 
rigorous  and  more  simple,  of  vector  densities  or  fiuxes,  which  is 
implied  by  the  field  equations. 

Returning  now  to  the  equation  (2.11)  we  see  readily  that 

the  ratio  of  the  first  to  the  second  term  on  the  right-hand  side 

is  again  / R (.  This  is  completely  analogous  to  conventional 
t 21  ) 

hydrodynamics  where  the  ratio  of  the  dynamical  to  the  frictional 

terms  is  There  is  this  difference  in  practice,  namely,  that 

we  can  on  occasion  realize  fairly  high  values  of  R in  the  laboratory, 

whereas  with  conventional  materials  R remains 

m 
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sm2.ll;  hence  while  we  are  familiar  with  turbulence,  we  are  not  well 

acquainted  with  the  type  of  phenomena  characteristic  of  cosmic 

hydromagnetism.  Noting  that  (2.11)  has  the  dimension  coB  \ we 

>-  j 

may  write  the  ratio  of  the  two  terms  on  the  right 

{ Rm  j = { V/Xa5dj  = ("v^d ) (2.16) 

where  G>d  is  a characteristic  frequency  of  the  free  decay  of  the 

field  in  tho  absence  of  motion  and  q is  characteristic  of  the 

v 

fluid  motions.  We  see  here  the  physical  basis  of  a dynamo  theory: 
for  sufficiently  large  the  magnetic  field  can  be  deformed 
(amplified)  by  the  fluid  motion  before  it  has  had  time  to  decay. 

It  remains  for  the  dynamo  theory  to  show  that  the  deformation  can 
occur  in  a sufficiently  ordered  fashion  so  that  a mean  magnetic 
field  can  survive. 

Next,  consider  the  mechanical  motion  of  the  fluid,  from 
conventional  electromagnetic  theory  we  have  for  the  density  of  the 
ponderomotive  force  of  the  magnetic  field 

_F  ^ _J  xJ3  = p,  X(  \7  x &_)  x B_  (2.17) 

The  corresponding  forces  of  the  electrostatic  field  are 
small  by  the  same  arguments  as  before.  The  force  (2.17)  will  appear 
in  the  Stokes-Navier  equations  for  the  fluid  motion.  Using  the 
well-known  vector  identity 

(V  xBj  x£  = (£•  V)B  - |\7(Bf)  (2,18) 

we  shall  write  these  equations  in  the  form 

ov/3t  + ( v . V7 ) v = -Vy  + (f  p)”1^  • V )B_  + V P’2v  (2.19) 
where  we  have  set 

= p + U + (2n)”1Bf  (2.20) 
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U being  the  gravitational  potential.  The  equations  (2.11)  and 
(2.19)  together  with  ouch  subsidiary  conditions  as  an  equation  of 
state,  determine  the  dynamics  of  a cosmic  fluid  in  which  magnetic 
fields  are  present. 

In  the  case  of  i n c © 9 sibility,  \/  • = 0,  this  system 

admits  of  a remarkable  symme tr iz a t ion  which  puts  in  evidence  the 

analogous  roles  of  and  v_  in  the  theory.  Letting 

1 I 

2 = v + (f  , J^  = v - 

2v,  = V + V , 2vo  a V - V 

j.  m e m 

we  obtain,  on  adding  or  subtracting  (2.11)  and  (2.19) 
sp/at  + (q  • \7)P  = +^72(v-,P  + v„q) 

ffA.  — * t -L  {.* /*w 

_ — P (2.21) 

ac^/at  + (^  • \7)q  = - + V (V;l9_+  Vgjp) 

where  now  (2.20)  becomes 

V'  » ^-1(P  + U)  + (2.-  5.)2/8  (2.22) 

These  equations  were  derived  by  Lundquisb  (1952)  independently  of 
the  author  (1950a).  The  symmetry  of  these  equations  mi6;4b  be 
somewhat  misleading:  thus  we  notice  that  Vg  becomes  negative  when 

electromagnetic  dissipation  outweighs  frictional  dissipation,  a 
fact  that  has  no  analog  in  ordinary  hydrodynamics.  The  ratio 

V/Vm  = Rn/R  = V (2.23) 

may  be  estimated  from  elementary  kinetic  theory  for  an  ionized 
gas  such  as  hydrogen  (El3asser,  1954).  One  finds  the  numerical 
value 

per  v ~ 2*10“4a^ 

where  a Is  the  degree  of  Ionization,  / the  density  in  mjes  units. 


This  shows  that  electromagnetic  dissipation  outweighs  frictional 
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loss  in  the  Interior  of  the  stars  where  f is  moderately  large, 
whereas  in  the  rarefied  intragalac tic  medium,  and  also  in  regions 
such  as  the  solar  chromosphere  or  corona  the  electromagnetic  loss 

^ 1 r QwTitfr!*  r~-  f ft  "]  a p « Vstt  t?j  ftftVifjft”?  ft  a 1 ■»/  ^ S f5  9 1 i^TT  ^ 

We  next  turn  to  the  conservation  law3.  It  may  readily  be 
shown  that  neither  energy  nor  vorticity  is  conserved  for  the  fluid 
motion  alone.  The  war’:  done  by  the  fluid  is  just  the  negative  of 
the  work  done  by  the  ponderomotive  force  (2,17).  The  power  del- 
ivered per  unit  volume  is  thus 

-v  • 7 = n-J"v  x B • (S7  x B)  (2.24) 


We  may  of  course  obtain  the  same  expression  from  the  magnetic 
field  equations.  On  scalar  multiplication  of  (2.11)  with  B and 
transformation  of  the  first  term  on  the  right  by  a well-known 
vector  identity  we  have,  neglecting  dissipation 

( 2Ji)"1SB2/3t  = • [( vxB)xB  + ^(vxB)  • (2.25) 


On  integrating  over  a volume,  the  first  term  may  be  cor. 7erted  into 
a surface  integral  and  can  be  made  to  vanish  if  the  surface  is 
extended  to  a region  where  v = 0 (this  term  represents  essentially 
the  Poynting  flux,  as  may  be  seen  by  substituting  (2.10)).  The 
last  term  of  (2.25)  is  identical  with  (2.24). 

The  curl  of  JP  does  not  in  general  vanish  and  so  there  is 
transfer  of  vorticity  between  the  fluid  and  the  field.  This  is 
most  conveniently  expressed  in  terms  of  the  Kelvin  circulation 
theorem  of  hydrodynamics.  Integrating  (2.19)  along  a closed 
contour  and  again  leaving  out  the  frictional  terra  we  obtain 


. dC  = ix-VdC  • (B  . S7)B 
= |i  V(\7  x B)  x B • dC_ 


(2.26) 


f 
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rt is  difficult  to  simplify  the  right-hand  3ide  further,  out  we 
shall  not  be  required  in  the  sequel  to  make  explicit  use  of  the 
theorem*  While  all  the  preceding  integral  theorems  fere  of  some 
interest,  the  one  most  significant  in  practice  is  the  conservation 
of  magnetic  flux  in  the  limit  of  small  dissipation,  given  by 
(2„I4).  We  shall  make  ample  U3e  of  it  later  in  the  application 
to  the  dynamo  theories*  for  further  analysis  and  illustrative 
examples  we  may  also  refer  the  reader  to  a comprehensive  treatment  by 
Lundquist  (1952). 

It  is  possible  to  derive  the  field  vectors  from  a vector 
potential,  though  the  relationships  are  slightly  different  from 
those  of  more  conventional  electrodynamic s . If  we  set.  as  usual 

B = A * E “ -3A/Bt  (2.27) 

the  first  equations  (2.1)  and  (2.2)  are  identically  fulfilled. 

The  second  of  (2.2)  or,  rather.  (2.8)  gives 

= v.  x (V  x (2.23) 

We  cannot,  however,  set  = 0,  since  the  divergence  of  the 

first  term  on  the  right  of  (2.28)  does  not  vanish  (see  Sec.  3). 

We  now  assume  R large  and  neglect  terms  of  order  R “1.  The 
m m 

last  term  of  (2.28)  being  of  this  order,  we  find 

9/3t(  V • £)  = V • [ v,  x (V  x A.)]  (2.29) 

By  virtue  of  (2.10)  which  also  holds  apart  from  terms  of  order 
R^  (2.29)  is  identically  fulfilled*  Finally,  is  deter- 
mined from  the  second  of  (2.1).  Thus  the  assumptions  (2.27)  are 
justified  in  the  approximation  in  which  electromagnetic  dissi- 
pation may  be  neg]e  cted. 

In  most  applications  it  is  possible  to  ignore  the 
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longitudinal  ( irrotational ) part  of  Ej  and  hence  of  A,  alto- 
gether, since  by  the  first  cf  (2.1)  it  does  not  give  rise  to  a 
magnetic  field.  If  a method  i3  given  whereby  A may  be  split 
consist ently  into  its  t n nn s v 0 r s 0 ( div”^gd!ic9"*fr,30  ) i t b 

longitudinal  part  (e.g.  on  using  a system  of  normal  modes.  Sec.  4) 
we  may  re-introduce  the  condition  SP' • A.=  0 to  supplement  (2,27). 
Since,  hov/ever,  the  divergence  of  the  first  term  on  the  right  of 
(2.28)  does  not  in  general  vanish  (see  e^u.  3.5}  wo  must  then 
supplement  (2.28)  by  the  condition  that  only  its  cransverse  com- 
ponent will  be  taken  into  account.  Since  the  decomposition  of  a 
vector  into  its  transverse  and  longitudinal  parts  is  linear,  this 
can  often  be  done  with  comparative  ease. 

In  the  applications  to  astrophys leal  pro-  ~.L  j.~~ 

blems  both  & and  R are  numerically  large.  This  means  not  only 
that  we  have  turbulence  but  also  that  there  will  be  an  entire 
hierarchy  of  eddies,  the  larger  eddies  feeding  energy  into  the 
smaller  one3,  according  to  the  usual  turbulence  theory.  The 
largest  eddies  correspond  to  the  largest  values  of  R possible, 
the  smallest  eddies  correspond  to  either  R 1,  or  Rm'-*-'  1 , 
depending  on  whichever  one  of  these  two  numbers  is  the  larger. 

If  there  is  intense  transfer  of  energy  between  the  mechanical  and  el- 
ectransgnetic  degrees  of  freedom,  at  least  among  the  smaller  eddies 
(Batchelor,  1950)  the  cutoff  of  the  turbulence  spectrum  must  be 
determined  by  whatever  mechanism  of  molecular  dissipation  is  the 
more  effective,  mechanical  friction  or  Joule's  heat. 

In  a turbulent  fluid  the  transport  of  physical  properties 
such  as  heat  content,  momentum,  vorticity,  and  so  on,  is  deter- 


mined by  the  corresponding  molecular  coefficients j the  same 


! 
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applies  to  thermal  dissipation.  Tims  we  siis.ll  be  forced  to  x op]  sc s 
Vin  an  eddy  magnetic  viscosity,  v * say,  or  else  c"  by  an  eddy 
electric  conductivity,  flrT  . It  is  well  known,  of  course,  that 
these  quantities  are  actually  functions  of  the  turbulent  -’ts  ts„ 

To  exhibit  more  fomally  the  magnetic  eddy  stresses  and 
the  magnetic  eddy  diffusivity,  we  remember  that  purely  mechanical 
stresses  can  be  expressed  as  the  divergence  of  a stress  tensor 
\ 3sg  for  instance  Somme rf eld,  1950).  We  shall  assume  in  the  re- 
mainder of  this  section  that  the  fluid  is  incompressible.  Wo 
shall  use  ter.sor  notation  in  cartesian  coordinates.  The  mole- 
cular viscous  stresses  are  the  divergence  of  the  tensor, 
v(dv,/'3x,_  + Bv,  /dx.).  Again.  ir.i  a turbulent  medium  we  have  the 

-x-  n.  XL  X 

Reynolds  stresses  which  are  derived  as  follows:  The  Euler 

equations  of  the  fluid  are 


li'rr 

W V 


i _ 


(2.30) 


Now  on  account  of  incompressibility  we  have 
9v\ 


ssr  vk 


+ V 
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We  now  set 


v,  = v,°  + v.^  (2.311 

x i i x 

where  v°  refers  to  the  "smooth”  and  v*  to  the  "turbulent”  com- 
ponent of  the  velocity.  We  choose  the  decomposition  (2.31)  so 
that^  on  averaging  over  the  irregular  motion^  v^4  = 0,  hence 
^7  = v^°,  and 


V V = 

ik 


+ vi  vk 


ij-  o c 

= Vi  Vk 


- S 


ik 


ik 


where  S 


is  the  Reynolds  stress  tensor.  If  we  insert  (2.31) 
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into  (2.30)  and  average  we  obtain 


5S 


ik 


Now  the  ponderomotive  force  (2.17)  may  be  expressed  as 
the  divergence  of  a Maxwellian  stress  tensor  (see  for  instance 
Stratton,  1941).  Referring  to  unit  mass,  this  relation  is 
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(2.32) 


where  is  the  usual  Kronecker  symbol.  We  now  split  into 

B1  = B1°  + B*  (2.33) 

where  the  same  conditions  for  the  averages  hold  as  for  (2.31). 

The  equations  of  motion  become 
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sf^lk  + T^)  (2.34) 

where  in  place  of  (2.20)  we  have  now  f^'=  p + U,  the  term  with 
2 

B having  been  absorbed  into  the  stresses.  We  see  that  the  mech- 
anical stresses  produced  by  the  tnr>hnl«nt  component  cf  the  magnetic 
field  are  exactly  analogous  to  the  purely  mechanical  stresses  of 
turbulence. 

A similar  transformation  may  '9  effected  for  the  magnetic 
field  equations  (2.11).  We  note  in  the  first  place  that 

(.Vx  "21  slr(viBk  ~ VV  (2.35) 

^ J K K 

The  parenthesis  on  the  right  represents  an  antisymme trical 

tensor  (it  being  well  known  that  any  vectorial  product  may  be  written 
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as  such  a tensor).  This  feature  distinguishes 
(2.11)  essentially  from  the  equations  of  motion 
tensors,  mechanical  as  well  as  electromagnetic. 
We  introduce  now  the  decompositions  (2.31)  and 


riEk  - Vi  = v 
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i ak  “ vk  Bi  + 1 


the  field  equations 
in  which  the  stress 
are  symmetrical. 
(2.33)  and  1st 


ik 


where  I is  the  antisymmetric al  turbulent  induction  tensor.  Now 


(2.11)  becomes,  on  omitting  the  term  in  v 
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(2.36) 


where  the  last  term  now  describes  the  turbulent  diffusion  of  the 
magnetic  field.  It  must  be  emphasized  that  in  spite  of  the  super- 
ficial similarity  to  the  last  term  in  (2.34)  there  is  the  funda- 
mental difference  that-  the  induction  tensor  is  antisymmetrical. 

In  any  one  situation  we  can  reverse  the  sense  of  energy  transfer 
between  fluid  and  field  by  merely  reversing  the  direction  cf  jr,  as 
is  apparent  from  (2.24).  Similarly,  the  sign  of  the  components  of 
I depends  on  correlations  between  the  components  of  v and  of  B 
and  the  sign  of  these  correlations  may  be  reversed  in  the  same  way. 
Whereas  the  stress  tensors  S and  T in  (2.11)  act  in  a way  quite 
analogous  to  molecular  viscous  stresses,  giving  rise  to  irrever- 
sible effects  only,  this  cannot  be  said  of  the  tensor  I?  the 
classical  proof  of  the  irre vers iblity  of  viscous  stresses  depends 
on  the  symmetry  of  the  stress  tensor.  For  stationary  isotropic 
turbulence  only  the  diagonal  elements  of  the  tensors  remain,  so 


that  I vanishes  in  this  case.  This  of  course  does  not  mean  that 
there  is  no  transfer,  but  that  the  transfers  in  the  opposite 
direction  balance.  We  can,  however,  have  a systematic  transfer. 
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and  of  either  sign,  if  the  turbulent  pattern  is  sufficiently  an- 
isotropic. ^his  question  is  closely  related  to  th at  of  eddies 
producing  feedback  in  dynamo  models  which  we  shall  take  up  in 

Sec.  7. 


We  may  inquire  into  the  comparative  magnitude  of  the 
tensors  S,  T,  I.  A number  of  authors  on  astrophysical  electro- 
dynamics have  assumed  that  equipartit ion  holds  for  the  energy,  at 
least;  in  erder  of  magnitude: 


(9..X7) 


1 ' 


In  this  case  it  is  readily  seen  from  (2.32)  that  {.T}  = t2j  = (Si- 
The  magnitude  of  I is  related  to  the  of  the  effective  eddies  in 
the  same  way  in  which  S and  T are  related  to  H.  As  Batchelor  (1950) 
points  out,  the  largest  eddies  are  driven  mechanically  and  their 
energy  is  degraded  before  they  can  create  an  equilibrium  magnetic 
field  of  their  own  dimensions*  and  of  magnitude  (2.37  )j  for  the 
eddies  of  smalls'*  dimensions  but  above  cutoff  we  may  assume  such 
equilibrium  to  prevail*  We  would  then  be  able  to  estimate  that  the 
ratio  of  the  eddy-stress  and  eddy-diffusion  terms  to  the  dynamical 
terras  is  of  comparable  order  in  (2,34)  and  (2.36).  There  is,  how- 
ever, no  ba3ic  need  tc  assume  that  the  equipartition  (2.37)  holds 
even  approximately  in  a rotating  system.  A turbulent  regime  is 
not  a statistical  equilibrium  but  a dissipative  process,  at  the 
best  a stationary  one.  For  an  equilibrium  the  existence  of  detailed 
balancing  is  a sufficient  (though  not  always  a necessary)  condition 
for  equipartition.  For  non-equilibria  the  deviations  from  the 
equilibrium  statistical  distribution  ( equipartition)  are  larger, 
the  more  the  system  deviates  from  detailed  baiarcing.  In  a rotating 
system  the  Coriolis  force  lacks  mechanical  reciprocity  and  hence 
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the  system  would  not  even  permit  of  detailed  balancing  in  a 
hypothetical  equilibrium.  Since  the  Coriolis  force  is  essential 
for  the  dynamo  mechanism  which  maintain'*  the  magnetic  fields,  it 
is  preferable  not  to  have  recourse  to  the  equipart ition  assumption 
(2.37)  but  to  tackle  the  problem  from  basic  dynamical  principles, 
which  is  what  we  shall  do  later  on. 


i 
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Electric  Fields,  Potentials. 


We  shall  now  investigate  more  closely  the  range  of  validity 
of  our  equations.  The  electromagnetic  field  equations  used  wore 
those  of  th.6  conventional  Maxwellian  electrodynamics  of  ponderable 
bodies.  There  are  certain  observed  cosmic  phenomena  which  cannot 
be  described  on  this  basis,  in  particular  radio  noise  and  the 
acceleration  of  cosmic  rays,  but  outside  of  these  there  exists  a 
strong  presumption  that  these  macroscopic  laws  should  hold  for 
the  slow  motions  ana  the  large  dimensions  of  t-b*  fluids  considered. 

As  before,  we  shall  assume  that  (3  is  numerically  small  and  Rm  num- 
erically large.  A question  which  immediately  pesos  itself  is  that 
of  the  magnitude  of  the  electrical  effects  associated  with  hydro- 
magnetic  phenomena,  end  this  will  now  be  considered. 

It  is  well  known  that  the  electromagnetic  field  equations 
of  ponderable  bodies  can  be  written  in  a relativist ically  invariant 
form  (Minkowski's  equations,  see  Von  Laue,  1921)  Here,  however,  we 
need  only  consider  the  terms  linear  in  P and  may  neglect  all  higher- 
order  terms  in  the  Lorentz  transformation.  The  kinematical  equations 
of  the  Lorentz  transformation  reduce  to  the  simple  form 

v ' - v - vQt,  t’  - t,  3/3t ' = 3/3  + v^  4 ^ (3.3) 


where  is  the  velocity  of  the  primed  system  with  respect  to  the 
unprimed  one  and  r the  radius  vector*  from  the  origin. 

In  the  same  approximation  the  field  vectors  transform  as 


Introducing  here;  the  assumption  that  Em  is  large  whicxi  may  be 
expressed  by  (2.10)  we  have 
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( TT  T?  /r£  l = f ) = ( Rft^  1 

1 ’o-'-  ; u"'o'~  j \~p  ; 

provided  v and  v are  cf  comparable  order.  Hence  the  transforms,  tic 


reduces  to 


S’  = E 


E 4-  V x B , B =B 

— ■ ~3  o ^w.  — 


(3.2) 
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vcA  EJ  = 8111,1  {<rE/j}  = (Rm| 

by  the  preceding  section,  we  see  that  is  small  of  order 

Now  Y i 3 so  exceedingly  small  that  this  product  is  small  for  all 


Pfinsrmohl  13  Vflbm.q  r>f*  R t Vionf'a 
-----  --  ~ m' 


JL'  - J*  > 1 


1 = h 


(3.3) 


the  second  equation  following  from  general  principles  of  relativity 
Furthermore,  the  conductivity,  or  , can  be  shown  to  be  a Lorentz 
invariant  from  general  thermodynamic al  considerations  ( fan  Lane,  3921 ) a 


wo  Oil  UU  l* 


that  all  the  equations  containing  alone 


will  not  change  under  the  transformations  considered,  nor  will  the 
ponderomotive  force  (2.24). 

We  next  inquire  into  the  space -charge , . Wo  must  have 

conservation  of  charge  which,  on  using  the  full  expression  (2.3), 
may  be  written 

“ it  = ^ • (£"  £ If)  = f 1 + B)  + V • (%  v) 

The  last  term  is  small  and  may  be  neglected,  and  we  are  left  with 
a differential  equation  for  7 f 

>5  + (<r  /c ) 7 = cr  f(t)  (3.4) 

where 

-f(t)  = • ( v x B ) = v • V7  xB-B  • V7  x v (3.5) 

^ AM  AM  <VV\ 

doss  not  in  general  vanish.  The  integral  of  (3.4)  is 


n (t)  = exp(  ~c  t/  £. ) / dt  f(t)exp(<r  t/c) 

' o 

Now  the  rate  of  charige  oi  f is  determined  by  the  frequencies,  to* 
of  the  fluid  motion,  and  c r /£  = Yo  <<  co.  On  letting  therefore 

f(t)  = f (o)  t tf’(o)  the  solution  becomes,  to  within  terms  of  the 
order  of  Y , 

>1  { t)  = £ f (o)  + £ tf  5(o) 

Its  meaning  is  apparent,*,  the  space  charge  is 

^ (v  x B)  (3.6) 

« /V*>» 

and  as  the  parenthesis  on  the  right  changes  with  time,  ^ follows 
this  change  quasistatically  to  within  terms  of  the  order  of  Y'  , 
that  is  synchronously  to  all  practical  purposes. 

While  thus  hydromagnetic  induction  does  in  general  produce 
a space  charge  in  the  conductor,  (2,1)  shows  that  the  associated, 
longitudinal,  part  of  the  field  does  not  affect  the  rate  of  change 
of  3.  With  a proper  choice  of  the  constants  of  integration  we  can 
say  that  the  longitudinal  component  of  J»,  whose  sources  are  ^ , is 
not  accompanied  by  a magnetic  field,  a familiar  result  of  conven- 
tional electrodynamics.  Now  we  have  seen  before  that  the  electrical 
energy  density  is  small  compared  to  the  magnetic  one,  and  similarly 
for  the  ponderomotive  forces,  ^or  this  reason  we  are  not,  in 
general,  interested  in  the  electric  fields  as  such  when  dealing 
with  problems  of  hydromagnetic  theory.  These  fields  become  of 
physical  interest  only  If  perchance  it  can  be  shown  that  somewhere 
in  the  universe  they  give  rise  to  the  acceleration  of  elementary 
particles,  a topic  which  Is  somewhat  beyond  the  confines  of  our 
present  subject.  Thus  it  is  altogether  legitimate  to  ignore  in 
the  sequel  the  longitudinal  component  of even  though  It  is  of 
the  same  order  of  magnitude  as  the  transverse  component.  If  we 


w 


introduce  « vector  potential  wo  m«7  correspondingly  set  v7  * A ~ 0. 
as  we  have  explained  in  Sec.  3,  although  this  divergence  is  by  no 
means  small  compared  to  x 

Cosmic  matter  is  practically  never  an  insulator  in  which 
static  charges  could  be  maintained  over  an  appreciable  length 
of  time  (the  ease  of  the  formation  of  atmospheric  thunderstorms 
by  electrostatic  effects  being  a notable  exception)!  hence  there 
will  be  no  static  charges  other  than  those  given  by  (3.6).  It  may 
be  convenient  to  introduce  static  charges  for  mathematical  reasons, 
e.g.,  at  the  boundary  of  a conductor  against  vacuum. 

It  is  seen,  then,  that  jg.  has  both  a divergence -free  and  an 
irrotational  component!  the  former  is  given  by  the  first  equation 
(2.1),  the  latter  can  be  expressed  by  (3.5)  and  the  relation 

(3.7) 

which  holds  to  within  terms  of  the  order  of  \ Both  com- 
ponents are  as  a rule  of  comparable  order  of  magnitude.  In  deal- 
ing with  problems  of  hydromagnetism,  especially  the  dynamo  theory, 
it  is  as  a rule  more  convenient  to  work  with  equations  that  con- 
tain J3.  al  one,  so  that  questions  concerning  the  electric  field 
become  irrelevant,  especially  since  the  ponderomotive  force  of 
the  electric  field  or  the  eleetris  stress  tensors  are  by  (2.7)  neg  = 
ligibly  small  as  compared  with  the  corresponding  magnetic  quanti- 
ties. In  the  hydromagnetic  phenomena  of  rarefied  gases  the  electric 
fields  generated  may  become  important  for  the  acceleration  of 
individual  particles  (cosmic  rays),  but  these  problems  are  not 
within  the  scone  of  the  present  review. 

There  is,  however,  one  point  that  might  be  touched  upon, 
namely  the  conception  advanced  by  several  authors,  that  the 
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conditions  for  the  acceleration  of  particles  from  thermal*  or  in 
any  event  small,  energies  to  higher  energies  are  particularly 
favorable  at  so-called  "neutral  points".  These  points  are  defined 
so  that  = 0 in  the  local  frame  of  reference,  that  is,  in  a system 
of  reference  in  which  \r_  = 0.  It  is  correctly  argued  that  in  such 
a system  the  particles  move  in  straight  lines  and  that  hence  the 
conditions  for  setting  up  an  electrical  discharge  are  much  more 
favorable  than  elsewhere  where  a discharge  is  quenched  by  the 
magnetic  field  (spiralling  of  the  particles  which  effectively 
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mean  field  &t  a neutral  point  and  E the  mean  field  in  the  fluid. 


we  see  from  (2.10)  that  Eq  is  small,  specifically 


(E0j  = (E/Rm|  " {>%} 


(3.8) 


Since  under  most  astrophysic al  conditions  Rm  is  quite  large,  the 
theory  of  discharges  at  neutral  points  (Dungey,  1953)  which  at 
first  sight  is  very  attractive,  needs  a thorough  revision  in  the 
light  of  this  last  result. 

Previous  to  the  development  of  hydromagnet ic  theory  with 
its  amplif icatory  mechanisms,  the  possibility  of  molecular  electro- 
motive forces  as  a cause  of  the  large-scale  electromagnetic  fields 
has  often  been  investigated.  Thermoelectric  potentials,  poten- 
tials caused  by  pressure  differences  along  a material  boundary,  and 
potentials  caused  by  differential  diffusion  of  negative  and  posi- 
tive ions  in  a density  gradient  are  typical  of  the  mechanisms 
invoked.  These  theories  (the  writer  pleads  guilty  to  having  once 
advanced  one)  have  all  had  this  in  commons  they  use  Ohm's  law, 

J - cr  F»  in  place  of  (2.8)  to  compute  the  currents.  Now  if  we 
replaoe  the  first  equation  (2.1)  by 


9B/et  = - \7x  (2.  + Ei) 

where  is  the  impressed  era!',  our  equation  (2.II)  becomes 

$B/dt  = -^7x  (v  x B ) + - V x E.  (3.9) 

In  a turbulent  medium  the  molecular  magnetic  diffusivity  must  be 
replaced  by  the  eddy  diffusivity  which  is  very  much  larger}  roughly 

{V)={Vm?  <3a0! 

By  (2.12)  this  means  that  the  conductivity  is  increased?,  or  the 
resistivity  decreased.  This  reduces  the  magnitude  of  the  fields 
generated,  both  by  the  first  and  by  the  last  term  on  the  right  of 
(3.9).  Since,  however,  the  existence  of  turbulence  requires  that 
there  are  large-scale  motions  on  which  the  turbulence  feeds,  the 
first  term  on  the  right  gives  rise  to  large-scale  hydromagnetic 
Induction  which  is  reduced  but  not  in  general  wiped  out  by  the 
turbulence,  as  we  shall  see  In  detail  later.  The  effect  of  im- 
pressed emf’s  is  different;  the  stationary  current  and  the 
stationary  magnetic  field  are  established  mors  rapidly  when  vm 
Is  larger,  but  the  field  produced  by  the  saturation  current  Is 
reduced  In  magnitude.  Prom  the  last  two  terms  of  (3.9)  the 
field  corresponding  to  the  stationary  state  is  of  the  order 

(BiJ  = (XEi./vm}  (3.11) 

and  according  to  (3.10)  this  field  is  smaller  by  a factor 
under  turbulent  conditions.  The  theories  referred  to  above  have 
had  great  difficulties  In  accounting  for  an  of  sufficient  mag- 
nitude to  explain  the  observed  fields,  anu  these  difficulties  are 
multiplied  if  the  turbulent  resistivity  is  introduced.  Most  mag- 
netic fields  of  large  dimensions  must  be  explained  by  dynamo  action 
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that  is  by  hydromagnatic  processes  in  which  the  field  has  been 
regenerated  many  times  over.  There  is  empirical  evidence  for  this 
in  the  fact  that  almost  all,  if  not  all,  stellar  magnetic  fields 
are  time-dependent  (Babcock  and  Cowling,  1953)*  Theoretically, 
we  must  say  that  the  long  decay  periods  for  magnetic  fields  in 
stars,  comparable  to  the  age  of  the  universe,  which  are  computed 
in  the  basis  of  the  molecular  vm  (Cowling,  1945)  must  be  reduced 
by  a factor  which,  in  stars,  is  likely  to  amount  to  a high 

power  of  ten.  Clearly,  for  a dynamo  theory  the  wav  in  which  the 
magnetic  field  originated  at  the  occasion  of  the  fir-  st  amp  13.fi- 
catory  processes  becomes  rather  irrelevant  as  compared  to  the 
mechanism  by  which  fields  can  be  regenerated  and  maintained. 


-34- 


4,  Free  eriodlc  Modes. 

The  integration  of  the  basic  equation  (2.11)  is  difficult 
even  with  a relatively  simple  geometry.  We  shall  proceed  to 
integrate  two  special  cases;  the  case  v = 0,  representing  free 
decay  of  a current  system,  essentially  in  a solid  conductor,  and 
the  case  vm  - 0 representing  pure  hydromagnet ic  induction.  In 
this  section  we  treat  free  decay?  from  (2.11)  and  (2.12)  we 
now  specialize  to. 

V2^  - ^ cr 33/at  = 0 (4.1} 

To  solve  this  equation  v/e  use  the  method  of  normal  modes  adapted 
to  the  aperiodic  case  (Elsasser,  1946/7)  which  differs  in  a number 
of  particulars  from  the  case  of  electromagnetic  oscillations 
(Stratton,  1941).  We  set 

JB(r,t)  = J^(r)exp(-y^t)  (4.2) 

and 

-''Y  = k^/ud*  = k^v 
' * m 

so  that  (4.1)  becomes 

= 0 (4.4) 

both  -/V  and  k being  assumed  real.  We  have  of  course  also 

+ k2E  =0  (4.o) 

From  (4.3)  the  time  of  free  decay  is  of  the  order 

<4-6> 

increasing  with  the  square  of  the  linear  dimensions,  v is  of 
the  general  order  of  unity  (inks)  for  metallic  conductors,  per- 
mitting a ready  estimate  of  the  order  of  the  decay  times. 
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tiroes  are  fictitious,  as  we  P ^ in  ted  out  before,  since  actually 

rj  the  magnetic  eddy  viscosity- 
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The  conventional  scalar  wavo  equation  can  be  solved  by 
separation  of  variables  in  11  different  systems  of  coordinates 
(Stratton,  1941),.  Prom  these  solutions  one  obtains  seta  of  ortho- 
gonal modes  by  imposing  suitable  boundary  conditions,  for  reasons 
too  complex  to  be  stated  briefly  the  method  cannot  readily  be 
extended  to  the  vectorial  wave  equation}  both  the  separation  of 
variables  and  the  boundary  conditions  present  special  mathematical 
difficulties , systoms  of  modes  can  be  found  for  three  types  of 
geometry*  plane  waves,  cylindrical  waves  and  spherical  waves.  The 
theory  has  grown  up  in  a more  or  less  ad  hoc  fashion,  though 
Stratton  attempts  some  systematization,  following  e arlier  work 
by  Mie,  Debye  and  Hansen.  We  have  tried  to  develop  the  method 
used  to  solve  the  vector  wave  equation  and  the  formalism  required 
for  boundary  conditions  and  orthogonality  in  a fairly  general 
fashion,  so  as  to  exhibit  both  the  scope  and  the  limitations 
of  the  method.  The  ensuing  formulas  may  be  used  in  all  the  known 
cases  of  solutions,  We  then  proceed  to  cons  true t^  for  spherical 
boundary  conditions,  the  actual  solutions  ana  to  discuss  in  detail 
orthogonality  and  normalization. 


Vector  Wave  Formalism, 


Tn  order  to  construct  solutions  of  the  vector  wavs 
equation  we  start  from  the  scalar  wave  equation 

+ k2^  = 0 ( 4 07 ) 

and  shall  derive  from  every  solution  of  this  equation  three 
associated  solutions  of  the  vector  wave  equation  (4,4)}  one  of 
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s of  the  form 
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^7  . 
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(4.8) 


where  i3  any  suitable  solution  of  (4.7  )>  these  are  the  long- 
itudinal waves  of  vanishing  curl-  A second  solution  is  of  the  form 

T = V x u^Z  (4.9  ) 

where  2^  is  a vector  field  to  be  specified  later.  A third  solution, 
finally,  is 

kgL=Vxx.  (4.10) 

The  vector  fields  (4.9)  and  (4.10)  are  transverse,  of  vanishing 
divergence.  There  seems  to  be  no  simple  way  other  than  (4.S) 
and  (4.10)  to  derive  a pair  of  linearly  independent  transverse 
vectors  from  a scalar. 

We  car,  rewrite  (4.9)  as 

T = xZf  yVxZ  (4.11) 

On  subs titu ting  this  into  (4.10)  we  find  after  some  straight- 
forward calculations,  using  (4.7), 

kj§L  ~ 7(Vr*  £) 

* Vv^(v  • z.)  - 2 * V)Z.+  x\7xj,  (4.12) 

Again,  since  and  a*’®  assumed  to  obey  the  wave  equation  (4.4) 
which  may  be  written 

V7  x V X 3 - k23  = u , v7’  x V X T - k2T  = 0 (4.13) 

we  have  from  (4.10)  and  the  second  of  (4.13) 

X7  x ^ = kT  (4.14) 

Equations  (4.10)  and  (4.14)  exhibit  a characteristic  symmetry 


between  the  vectors  S and  T 


Next,  from  (4.9)  and  (4.14) 


On  che  other  hand,  the  curl  of  (4.12)  agrees  with  (4.15)  only 
if  the  last  three  terms  In  (4.12)  are  the  gradient  of  a scalar. 

We  shall  consider  ^ as  a quantity  determined  by  the  particular 
boundary  conditions,  hence  related  to  the  system  of  coordinates 
chosen,  whereas  is  any  one  of  a set  of  orthogonal  modes.  For 
the  present  purpose  we  may  therefore  assume  to  be  effectively 

arbitrary,  and  7^  and  as  independent  of  each  other.  It 
follows  that  the  last  term  in  (4.12)  must  be  the  gradient  of  a 
scalar  in  itself  since  it  contains  , whereas  the  two  other 
term 3 contain  the  derivatives  of  <4^  . These  two  terms  both  in- 
volve ^ , they  have  arbitrary  but  in  general  different,  dir- 

ections; hence  each  of  them  must  again  be  the  gradient  of  a scalar 
separately.  From  the  first  and  the  third  of  the  throe  terms  there 
follows  7 • Z - const-,  and  \7  x JZ  = const.,  respectively)  hence 
Z is  a linear  function  of  the  cartesian  coordinates.  Combining 
this  with  the  requirement  that  the  middle  term,  ( v7 ,'/v  * ^ ):§_ 
must  be  a gradient  we  find  two  solutions:  One  is  clearly  Zl - 

const.;  we  may  for  instance  take  2^  to  be  a unit  vector  in  the 
z-direction.  The  other  solution  is  Z = r,  where  r is  the 
radius  vector  from  the  origin.  This  latter  choice  gives 
(^7 4^  • ^7  )r  = The  former  choice  is  used  with  plane  and 

cylindrical  waves,  the  cylinder  axis  being  in  the  z-direction; 
the  latter  is  appropriate  for  spherical  waves.  We  then  have  from 
the  preceding  formulas,  for  constant  .Z 

T=Vy  X&;  5.  = k'f£  + k^VO^/az)  (4.16) 

and  f or  J5  - in  polar  coordinates 
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X ~ v y*  x.  + k"“^7(  &rv^/or)  (4.17) 

We  may  summarise  the  last  two  formulas  by 

£=  + k"x^(y  )*  (4.18) 

We  shall  now  write  down  formulas  from  which  the  orthogonal  it 
relations  of  these  vector  modes  may  be  derived  for  appropriate 
boundary  conditions.  Individual  scalar  modes  which  satisfy  (4.5) 
will  be  distinguished  by  indices,  • . .j  the  vector  modes 

derived  from  by  (4.18)  and  (4.8)  will  be  designated  by 
S^,  Tx,  U^.  We  new  have 

/%  • T2dV  = * V^2  • ^ dV  = /V  x (sf2  Vvf/2)  • & dV 

( 4. 1.9 ) 

= / x zJdV  = /^V^2  x„£  *_n  dcr 

v/here  n_  is  the  unit  normal  to  the  surface,  directed  outwards. 

Next  we  find,  using  Green's  theorem 

•Oil  • s3dv  = ' 2 4V  t k2-x  /Vvi  • V(^2)'dv 

(4.20) 

= Wa^Vi  • 2 dv  + k2'lkl2/‘tl'+/2dV  • k2'bV2),OVi/9n)d0- 
and  finally 

/Si  • ~VV  “ fcg'Vv^)'  x Vyj  dv  (4>21) 

= VVV-ftvf-aJVyi  x zlav  = * z,  • n.d<r 

Some  further  formulas  are  required  for  the  orthogonality 
relations  among  vectors  of  the  same  type.  For  the  U-vector s we 
can  transform  the  integral  over  by  means  of  the 

conventional  Green's  theorem  as  in  (4.20).  For  the  T-vectora  we 
find  from  (4.18) 

/J1  • JgdV  = AV^1rV^’2)Z2dV  - /(Z  •\7f'1)(Z  • Vv^gidV  (4.22) 
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Further  ugefu  1 r « 1 a l.  lulls  follow  from  the  vectorial  equivalent  of 
Green's  formulas  (Stratton,  1941).  For  two  arbitrary  vectors, 
provided  \7  * J3  = 0. 

/[{  \7  x A)‘{\7  x J3 ) *■  = fjk.  x ( V x Ji ) • ndCT  ( 4.23 ) 

which  may  readily  be  proved  on  converting  the  right-hand  side  into 
a volume  integral  by  Gauss'  theorem,  ^he  equivalent  of  Green's 
second  formula  is 


/[*SV2£’£'V2£ ]-  = f[k  * (V  x B)  - J x ( V x^A)  j*  dn«r(4.24 ) 

provided  K?  ' A.  ~ ^ * JB  - 0.  From  (4.23)  we  obtain  an  expression 
relating  the  orthogonality  of  the  S-vectors  to  that  of  the  T-vectors 


k-.k^S,  • S0dV  = v'T,  * TodV  + k^S,  x T„  * ndcr-  (4.25) 

± N*  -L  N*  O CZ  _L  fWK  6 G>,  w\  _L  «M  C* 

The  surface  integral  on  the  right  can  be  further  reduced  by  means 
of  the  identity,  following  from  (4.18) 


./fc,  x T0  • n d<r  = /P0  • n x S,  dcr 

«*i  X A"uc!  WX<£  <*«  ft^X 


=:  VK1(  ^V2  X-S)*(JL\  x Z)d<r 


(4.26) 

kjV(VY2 


From  (4.24)  we  get  the  two  useful  relations 

(k®  - k|)Ax  • JgdV  - /[k^  :t  ig  - k,J2  X sj  • n dcr<4.27) 

- k|)/§1  • J2<JV  = /[kgSj  x t2  - k1§2  x • n art  4.2S> 

If  our  conducting  body  is  of  finite  size  (rectangular  box 
or  cylindric  prism  for  cartesian  coordinates,  circular  or  ellip- 
tic cylinder  for  cylindrical  coordinates,  sphere  for  polar 
coordinates)  we  must  join  our  solutions  for  the  inside  to  the 
solutions  of  Maxwell's  equations  for  free  space  at  the  outside. 

Assuming  that  6.  and  *x  have  the  same  values  inside  the  conductor 
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as  outside,  the  boundary  conditions  impose  continuity  of  all  field 
vectors,  except  that  there  may  be  a surface  charge,  "TT,  per  unit 
area  so  that  for  the  normal  components 


( E ) 

v n outside 


*"^n  ^ in3ide 


(4.29) 


In  outer  space  ’.ve  have  a field  wjfttich  adjusts  itself 
qunsistatically  to  the  fields  that  prevail  at  the  boundary  and 
that  are  the  result  of  the  validity  of  (4;1)  and  (2.1)  on  the 
inside.  Since  the  displacement  current  is  negligible,  we  have 
from  (2,1)  and  (2,2)  for  the  outside  field 

^72ML  = o , T72 E = 0,  V*JBl=0,  V • £=  0 (4.30) 


It  is  not,  in  general,  possible  to  choose  boundary  con- 
ditions for  the  solutions  of  the  vector  wave  equation  such  that 
a fully  orthogonal  system  of  vectors  J5,JT,  results.  The  boundary 
conditions  on  are  of  course  determined  by  the  boundary  con- 

ditions of  the  electromagnetic  field  quantities;  it  then  appears 
that  "almost”  all  the  modes  are  mutually  orthogonal,  but  that  for 
certain  pairs  of  vectors  the  orthogonality  fails  (Stratton,  1941). 
Fortunately,  this  failure  is  sufficiently  limited  so  that  it  does 
not  seriously  hamper  the  use  of  the  formalism  for  the  solution  of 
our  physical  problems.  There  is  one  case  where  full  orthogonality 
obviously  obtains:  for  plane  waves  in  a rectangular  box  with 

cyclic  boundary  conditions  at  all  faces. 


Modes  of  the  Sphere. 

We  shall  be  principally  interested  in  spherical  conductors. 
We  assume  uniform  conductivity  and  let  R be  the  radius  of  the 
sphere.  We  shall  assume  the  outer  space  to  be  vacuum.  Consider 
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first  the  inside.  W'e  have  for  the  solutions  of  (4.7), 
from  an  arbitrary  constant. 


'r  - *> 


where  j is  a spherical  Bessel  function. 


J_(x>  = (it/2x)1/2J 


-n 


1+1/2 


(x) 


apart 


(4,31) 


(4.3a) 


and  where  the  k are  determined  by.  the  boundary  conditions* 
ns 

•The  are  the  conventional  surface  harmoncis.  we  now  writs 
down  the  expressions  for  the  vector  modes  in  polar  coordinates  in 
terms  of  vXx  « They  are  given  by  (4.8)  and  by  (4.17)  which  yield 

Jr  » d‘-p/ar  , = r_1e^/c^ 

U = (r  sin  i^-)’”~c)'^/df 


(4.33) 


This  type  of  mode  13  purely  longitudinal.  Next  we  have 

This  type  of  mode  will  be  designated  as  toroidal.  Finally 
JSr  = kr  + k-1a2(r  v^)/9r2  = n(n+l)  (kr)”1i^JN' 


(4.34) 


3 = (kr  )”^82(r  y)/drdifl- 


} (4.35) 


J5  = (kr  sin  iP-) ’^c*2(r  v^J/Qro^ 


i 

J 


This  type  of  mode  will  be  designated  as  pololdal. 

Ffv»  cm  mif.ATi  flrio  r>  c*  f 1 a!  H Artliof  I r-.nfl  nonnoo  -r.  n 4*.  Vio 

* - £ ~ - - 

vectorial  Laplace  equations  (4.30)i  the  corresponding  generating 

scalar  is 


V = c r"n-1y®(^-,  <p) 


i A *Z 
\ Q 


\ 

/ 
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wher-e  C trill  be  chosen  so  a3  to  assure  continuity  with  the  inside 
solution  by  virtue  of  the  electromagnetic  boundary  conditions. 

The  solution  defined  by  (4.8)  or  (4,33)  continues  to 
hold  in  the  limit  k = 0.  In  the  case  of  the  transverse  modes, 
however,  we  must  redefine  our  expressions  sc  that  (4.35)  remains 
finite  for  k = 0.  This  is  easily  done  by  considering  the  vectors 
kS,  and  k£  in  pl&ce  of  the  above,  ^hen  kT,  vanishes  as  k goes  to 
zero  and  one  readily  verifies  that  kS^  = -nU>  this,  furthermore, 
agrees  with  (4,14),  namely  £ = V'  x (,§/k)  ~ 0.  Hence  on  the 
outside  there  exists  only  one  solution  which  is  essentially 
longitudinal,  but  may  if  desired  be  expressed  in  terms  of  a 
pololdal  vector  field  with  the  generating  function  (4.36). 

We  shall  now  construct  the  free,  apericdical ly  damped 
electromagnetic  modes.  Again,  will  d esignate  field  vectors 

from  which  the  time  factor  has  been  split  off,  by  (4.2).  The 
longitudinal  modes  (4.33)  are  readily  constructed,  but  they  are 
purely  electrostatic  and  the  corresponding  magnetic  field  vanishes. 
As  indicated  in  Sec,  3 they  are  only  of  subordinate  interest 
for  the  dynamo  problem j their  formalism  is  moreover  quite 
s t r aight forward • 

For  the  transverse  modes  we  have  from  the  field  equations 
together  with  (4.2)  and  (4.3) 

V x 3^  = k(iiC-jS/k) ; , \7  x (uo’.g/k)  = kB^  (4.37) 

If  we  introduce  a vector  potential  it  is  related  to  S by 

A = E/£A.  = ii(T E Az*' 

' 

The  equations  (4.37)  are  identical  in  form  with  (4,10)  and 
(4,14).  Moreover,  the  vector  wave  equation  (4.13)  is  a 
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consequence  of  these  r3l3.olcn3«  * ne r o arc  two  b y p e s of  transverse 
modest  the  toroidal  magnetic  modes  where  J?  is  of  typo  jT.  and  the 
poloidal  magnetic  modes  where  J*  is  of  type ^3, 

We  first  consider  the  toroidal  magnetic  modes.  To  obtain 
an  idea  of  their  geometry  we  note  from  (4.34)  that  for  rotational 
symmetry  (zonal  harmonics)  there  is  only  a cp-component  j hence  the 
magnetic  lines  of  force  coincide  with  the  circles  of  latitude. 

For  the  toroidal  dipole  mode  has  the  same  sign  throughout, 
whereas  for  the  toroidal  quadrupole  mode  J3  changes  sign  at  the 
equator-. 

We  have  seen  that  the  .T-vectors  vanish  identically  in 
outer  space  where  cr-  =•  0.  Hence  the  boundary  condition  is 
J3  = 0 at  the  surface  of  the  conducting  sphere,  r = R,  From 
(4.34)  and  (4.31)  this  gives  the  characteristic  equation 

VknsR>  " 0 <4*S9> 

The  k form  a twofold  sequence,  depending  on  the  "quantum  numbers" 
n and  s,  the  latter  numbering  the  successive  zeros  of  The 

electric  field  corresponding  to  this  mode  is  poloidal)  by  (4.35) 
its  normal  component  vanishes  at  r = R,  but  the  tangential  com- 
ponents do  not  in  general  vanish.  We  can  fulfill  the  boundary 
conditions  by  combining  an  external  multipole  field  with  a surface 
charge  defined  by  (4.29  )„  (It  3houlu  be  pointed  ou  t that  the 
transverse  vector  modes  start  with  n = lj  the  solutions 
corresponding  to  n = 0 vanish.) 

We  xie A.t  consider  the  poloidal  magnetic  modes.  For  zonal 
harmonics,  that  is  rotational  symmetry,  the  magnetic  field  lines 
are  confined  to  the  meridional  planes.  T0  fulfill  the  boundary 
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condition  for  3^  for  the  general  mode  of  this  type  we  must  have 
continuity  of  all  three  components  at  the  surface  of  the  conductor. 

It  is  convenient  to  express  the  external  multipole  field  in  terms 

1KL 

ofA  poloidal  vectors  (4.35)  rather  than  by  the  longitudinal 
U.-vectors?  we  have  shown  above  that  this  is  legitimate.  Now  for 
this  external  field  whose  generating  scalar  is  (4.36)  the  relation 
of  the  tangential  to  the  normal  components  of  j§L.  is  given  by 
9(rv|/)/dr  = -n  ^ . If  there  is  to  be  continuity  of  all  components, 
the  internal  must  obey  the  sajr.s  condition  at  the  surface, 
r = R.  that  is 


3(rJ  ) 

— 55—  + nJn 


**  0 


(4.40) 


= R 

This  is  the  characteristic  aquation  for  these  modes  which  may 
readily  be  transformed  into 

.1..  , (kP_R)  = 0 (4.41) 

-xi-x  ns 

Moreover,  If  (4*36)  is  the  generating  scalar  for  the  external 
S-vectcrs  we  obtain 


C = p'n+-M  (1r  D) 
ns  n •'n  na  ' 


(4,42) 


The  electric  field  of  these  modes  Is  purely  toroidal.  Since  there 
can  be  no  toroidal  field  In  empty  outer  space,  the  boundary  con- 
ditions for  J*.  cannot  be  fulf  illedj  In  order  to  satisfy  them  we 
would  have  to  go  to  a higher-order  approximation.  Since  the 
aperiodic  modes  decay  extremely  slowly,  we  could  try  to  fulfill 
the  boundary  conditions  by  assuming  a minute  electrical  conduc- 
tivity in  outer  space  which  would  make  a toroidal  electric  field 
possible  and  would  also  correspond  fairly  closely  to  geophysical 
and  astrophysical  conditions.  In  an; 


event,  the  electric  field 
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In  again  negligibly  small. 

We  next  discuss  the  orthogonality  and  normalization  of 
the  modes.  From  (4.20)  and  (4.21)  one  readily  verifies  that  all 
T - vectors  are  orthogonal  to  all  J5 -vectors  over  the  interior  of 
the  sphere.  Similarly,  it  follows  from  (4.19)  that  all  Ty vectors 
are  orthogonal  to  all  Uyvectors,  The  J3-vectors  are,  however,  not 
orthogonal  to  the  JJ- vectors  (Stratton,  1941)  bnf  since  we  are  not 
interested  in  the  longitudinal  field  components  this  fact  will 
not  impede  our  calculations.  “Next,  consider  the  mutual  orthogon- 
ality and  normalization  6f  the  toroidal  modes.  introduce  the 

abbreviation  S?  J by 

V = r_1r(a/ar)  + r^V*  ( . (4.45) 

so  that  ^7 ' is  a gradient  vector  along  the  surface  of  the  unit 
sphere.  Applying  (4.43)  andl  (4.31)  to  (4.22)  we  find 

• T2dV  = /J1j2r2dr  . ^Yj.dCT  (4.44) 

To  evaluate  the  surfaco  integral,  consider  for  a moment  a potential 
function,  ^ = r11^.  Applying  Green’s  theorem  to  two  such  func- 
tions we  find,  on  using  (4.43), 

o . d'+'o 

/r'Wi  • V-  y 2dv  + / -gS.  dv  = /H^1  acr 

If  we  now  choose  as  the  volume  of  integration  ths  interior  of  the 
unit  sphere,  the  integrals  are  readily  evaluated  and  yield, 
provided  we  use  complex  harmonics,  = P^(cos 

X = / Y_  • T7':y_dcr  = o for  Y-  ^ Y*  (4.45) 

x id  x ' a 

where  the  asterisk  designates  the  conjugate  complex,  and 

x = / = nin+l)  /|  ^!2dcr 

_ 4rcn(  n+1)  (n-Ha)  1 
2n  + 1 " ( n=m ) J 


( 4.45 ) 
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In  order  to  prove  orthogonality  of  the  radial  functions 
for  fixed  n and  variable  x we  require  the  boundary  condition  for 
the  toroidal  modes  which  is  0 at  the  surface  of  the  sphere, 

by  (4.39),  With  this,  the  radial  integral  of  (4*44)  may  be 
evaluated  from  the  formula,  proved  in  the  theory  of  Bessel 
functions  (see  Jahnke-Emde , Tables  of  Functions) 

/ J^(icx)x2dx  » (x3/~)[ J^(kx)  - jr„i(kx) jn+1(kx)7  (4.47) 

which  now  reduces  to 


-<R3/2>W^sH)Jn+1(k>) 


(4.48) 


The  right-hand  side  is  of  course  essentially  positive  in  spite 
of  the  minus  sign  in  front. 


For-  the  puioiaai  modes  we  again  discuss  first  orthogonality 
with  respect  to  the  spherical  harmonics.  Using  (4.43)  in  (4.17) 
we  see  that  S = arY  + b ^7  *Y  where  a and  b are  functions  of  r, 
and  orthogonality  with  respect  tc  the  Y*s  follows  from  (4.45). 

To  establish  the  orthogonality  of  the  radial  functions  and  the 
normalization,  we  first  derive  from  (4.28),  (4.26)  and  (4.43) 
the  relation 


X 

TT 


R^kl"k2^ 


% 

it  J1 

k2  x 


5r 


k2 

*1 


3(1*;^) 

3r 


(4.49) 
r = R 


where  X is  the  surface  integral  (4.45-46).  If  we  multiply  the 
bracket  by  k,kg  and  substitute  the  boundary  condition  (4.40)  the 
right-hand  side  is  seen  to  vanish  for  n,  = provided  however, 
k^  f kgj  for  k^  = kg  the  formula  becomes  invalid.  Thus  the  pol- 
oidal  mod  63  are  orthogonal  uvoi-  the  interior  of  the  sphere.  Since 
for  given  n the  different  radial  modes  have  the  same  external  field 
apart  from  constant  factors,  the  modes  are  not  orthogonal  on 
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integration  over  all  space.  For  k,  — k„  we  use  (4.25)  and  (4.26) 


and  obtain 


A l jm  ! 2 A i _f  2 1 

/ . a ! av  - J i Tt  av !y 

& 


3(PJ„) 


n 


3r  . 


i • • ^ 

I V'M  cur 


r=B 


The  first  integral  on  the  right  is  given  by  (4.44)  but  to  evaluate 
it  by  (4.47)  we  must  now  use  the  boundary  condition  (4.41).  At 
the  same  time  we  can  simplify  the  bracket  on  the  right  by  im  ans 
of  the  form  (4.40)  of  this  boundary  condition.  The  result  is 

/ U/2dV  = |r3/2  + nR/knSl^(knaR)  / / ^’y/^cT  (4.50) 

inside  L J 

We  finally  compute  the  overlap  of  the  external  fields  for 
the  same  n but  different  values  of  s.  Since  most  of  the  preceding 
general  formulas  break  down  if  applied  to  the  outer  space,  the 
calculations  are  conveniently  carried  out  directly  from  (4.55)  on 
using  the  generating  scalar-  defined  by  (4.36)  and  (4.42).  The 
result  is 

/|  3pdV  = (nRA^kgU  1(k1R)j2(k2R)  /\\7'Y\2dCT  (4.51) 

outside 

For  k.  * k„  this  corresponds  exactly  to  the  second  term  in  the 
bracket  of  (4.50).  These  expressions . represent  of  course  twice 
the  magnetic  energy  of  such  a mode  on  the  inside  and  the  out3id6, 

respectively.  Hence  the  ratio  of  the  outside  to  the  inside  energy 
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is 


!outside/Einside  f1  + ^knsR^  /2nJ 


(4.52) 


Let  us  finally  obtain  an  estimate  of  the  decay  time  of 
these  modes  for  the  earth’s  metallic  core.  We  shall  see  later 
that  the  toroidal  dipole.,  modes  cannot  be  appreciably  excited, 
the  most  significant  modes  being  the  poloidal  dipoles  and  the  tor- 
oidal qucdrupoles.  For  the  former  we  have  j (kR)  = 0 and  for  the 
latter  = 0.  The  lowest  root  for  the  dipole  is  kR  - n,  the 

lowest  root  for  the  toroidal  quadrupcls  is  kR  ~ 5.8.  From  (4.2) 
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and  (4*3)  the  decay  time  .13  uer/k'5.  With  R = 3.5  • 10^  meters 

Q 

and  O'  = 10  mks,  one  tenth  of  the  conductivity  of  ordinary  iron, 
we  obtain  50,000  y®ars  and  14,000  years  respectively.  These 
figures  are  of  course  purely  nominal,  since  the  actual  decay  times 
are  determined  by  the  magnetic  eddy  diffusivity  and  are  no  doubt 
much  smaller,  perhaps  closer  in  order  of  magnitude  to  1000  year3, 
as  may  be  judged  from  certain  features  of  the  geomagnetic  secular 

IT  a *r»  -?  q ¥ •?  on  / T*-!  ~ *1  c»  C /"v  % 

^ cu  X S t — Uxl  \ a-  y j.«7uv/i 
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5.  Kinematics  of  Induct  Ion, 

Ws  now  consider  the  induction  equation  in  the  absence  or 
dissipative  losses;  by  (2.11) 

9B  / 9t  = S7  x ( v x B ) (5*1) 

AM 

or  its  equivalent  (2.15),  namely, 

d(5*  ”1B)/dt  = (^“1B.  • V)v  (5.2) 

and  the  corresponding  integral  theorem  (2.14)  which  is 

d/dt  / = 0 (5.5) 

In  terms  of  the  vector  potential  we  have  from  (2.28) 

aij/at  = v^x  (\7  x a)  (5.4) 

which  may  be  written  alternately 

-A/dt  = vx(\7xA)‘l*(v»V)A  (5.5) 

^ AA.  OMA  "VMA 

and  (5.3)  becomes 

d/dt  / A,  • dC.=  0 (5.5) 

It  is  at  once  apparent  that  the  induction  process  does 
not  involve  any  material  constants  of  the  medium.  The  relative 
rate  of  change  of  the  magnetic  field  is  of  the  order  £ v/A  j ; 
ahus  if  we  wait  for  a time  during  which  a fluid  particle  travels 
a distance  of  the  order  of  the  linear  dimensions  of  the  system, 
the  amplification  of  the  field  can  become  appreciable;  for  much 
longer  times  the  amplification  may  become  very  large  under  con- 
ditions otherwise  suitable  which  will  be  discussed  later. 
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It  has  been  found  by  Parker*  (1954)  that  the  conservation 
equations  given  above  can  be  integrated  with  respect  to  the  time. 
In  order*  to  avoid  some  of  the  mathematical  complexities  of  the 
full  tensor-analytical  treatment  we  shall  adopt  a mixed  formalism 
In  which  we  pass  from  vector  to  tensor  notation  as  required. 

We  integrate  (5.2)  over  a volume),  getting 

/ f (d/dt)(B/f)dV  = / (B  • V)z.dV  (5.7) 

Now  we  may  show  that 

J (B  * ^7)v  dV  = J v B do-  (5.8) 

~~  — — n 

To  prove  this,  consider  the  x-somponentj  we  have  the  identity 

JL  * ^vx  = ^ * (vx§>)  " vx^  * ^ 

and  the  last  term  vanishes.  On  applying  Gauss'  theorem  on  the 
right,  (5.8)  follows. 

Let  us  assume  that  the  volume  of  integration  Is  attached 
to  the  fluid  particles  and  moves  with  the  latter.  Then  Bnd^“ 

r\ 

Is  an  Invariant  by  (b.3)  and  may  be  written  (Bnd<r,)~  where  the 

superscript  o will  here  and  in  the  sequel  designate  the  values 

of  variables  at  an  initial  Instant,  t = 0.  -hus  we  shall  write 

/B° , v°,  r°  for  field,  velocity,  and  position  pertianing  to  a 

particle  at  time  t = 0-  This  is  essentially  the  Lagrangian  method 

of  hydrodynamics  where  the  variables  referring  to  time  t are 

considered  as  functions  of  the  variables  which  characterize  the 

same  fluid  particle  at  time  t = 0.  Thus  B.  v.  r are  functions  of 

B°,  v°,  r^  and  of  in  components  we  consider  B . , v.,  x4  as 
***  *“  X X 

functions  of  the  initial  values  B°,  v°,  and  of  the  time. 
simplify  the  fo:malism  we  shall  confine  ourselves  for  the  present 
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to  cartesian  coordinates.  We  writs  for  the  operator  with  the 

components  o/dx°. 

Now  since  B do-  = B°dC  0 by  (5.3).  we  see  from  (5.8)  that 
n n J 

/ (B^  * \7)v  dV  = J (_B°  • V7°)v  dV° 

where  on  the  right  v is  considered  a function  of  the  x°.  Since 
y d'v  = j°dV°  this  gives 

/ L(B  • V)  - (f/r°)(B°  • \7°)jv  civ  = o 

Here  we  can  equate  the  integrand  to  zero  since  the  volume  of 
integration  is  arbitrary.  If  we  substitute  in  (5.7)  we  find 

d/dt(B/f)  = [(BC//°)  • V°jv,  (5.9) 

This  integrates  to 

B/f  - B°/f°  - [iB7f0>  • V°]ir-r°) 

But  since  B°  = (B_°  • O'  ° )r°  . identiallly,  this  reduces  to 

S/f  = [(B0//0)  • V°]r  (5.10) 

which  is  the  desired  integral,  expressing  at  time  t in  terms  of 
E>°  and  the  kinematical  properties  of  the  fluid. 

We  next  discuss  the  conservation  theorem  for  the  vector 

* } 

potential.  Equation  (5.5)  is  closely  related  to  a well-known 
vector-analytical  identity.  If  we  let 

^7(v  . A)  = [\7<V  • +[V(v,-  *)]A=oonst.  (5.11) 

&3  find  readily  that  (5.5)  may  be  written 


d£/dt  = [V(v  • A )J 


v=const . 


(5.12) 


For  several  of  the  subsequent  formulas  of  this  section  the 
author  Is  indebted  to  Dr.  Willaim  L.  Bade. 
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If  we  substitute  (o.l2)  into  (5.11)  we  get 


dA/dt  = -A  x (V  x v)  - fA»V)v  + ^7'4-'  (5.13) 

where  v-P  = v • A may  usually  be  ignored  for  the  reasons  explained 
in  Sec.  3.  These  expressions  may  be  greatly  simplified  on  using 
tensor  notation.  By  viture  of  (5.11)  the  relations  (5.12)  and 
(5«13)  become 


r?  A / /3 


,»  f * ft  / V,  \ 

w'n’k/  w-f“i / 


-i:\OV9xp  + aY/a):i 


(5.14) 


To  obtain  the  time -integral  of  the  conservation  equation  it  is 
best  to  start  over  again  from  (5.6)  which  may  be  integrated 
directly  giving,  in  tensor  notation. 


But 


21  / (A,dci  - Ajdoj)  = 0 


(5.15) 


dc°  ~ HdcvOx5/B^) 

x k 

On  substituting  this  into  (5.15)  the  integrand  must  be  the 
gradient  of  a scalar? 

= XI  a£(  Bx^/dxn)  + (5.16) 

k 

which  is  the  desired  integral  for  A*  1 t may  be  compared  with  the 
corresponding  formula  (5.10)  for  _B,  which,  in  the  present  notation, 
reads 

B±/f  = ?EI(B°/P0)OxA/ax^)  (5.17) 

If  we  drop  the  gradient  term  ifn  (5.16),  we  obtain  a simple  relation 
on  forming  the  scalar  product  of  A~  and  JB/f , namely^ 
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i Tl  f & A ^ 

A • o/p  — A 


3°/f  ° , ( d/dt ) ( A * B ft)  = 0 


(5.18) 


So  far  we  have  confined  ourselves  to  cartesian  coordinates. 
For  reference  purposes  we  give  the  expression  of  the  preceding 
integral  theorems  in  cylindrical  and  in  spherical  polar,  coordinates. 
We  shall  omit  the  density/  ; to  reintroduce  it,  it  ’ s '•■erely  nec- 
essary to  replace  B by  B If  srtd  E;°  by  B °/f°.  The  calculations  are 

straightforward  differential  geometry  and  need  not  be  described. 

n 2 P 1 /2 

In  cylindrical  coordinates,  $ , 9,  z (where  5 = (x  + y'y'  is  of 
course  not  to  be  confounded  with  the  density)  the  field  transforms  as 

B = B?  + i-  14  p°  + *L  Bc 
' af°  * f°  S{p^  * sz°  ~ 


B - 


= 


4.  JL  B°  + f ^ n° 

5 3/°  3z°  55 


3z 

af° 


1 3z  ~o 


a«? 


ii  +■ 

c <p 


32 

3zl 


B. 


( 5.19) 


and  the  vector  potential  as 


.o  . Oo  3 

1 j 


IV  ^ 


3z°  o 
37"  Az 


A = i 8fl  A"  + £ A°  + i |5* 

% 7 V f \ f tut 


(5.20) 


. - 3f°  ac  4.  P°  4.  9Z°  A° 

Az  " 3T*  + y 82  % + 3z~  Az 


In  spherical  polar  coordinates,  r,  <p  we  rind  for  the  field 

1 3r  o A 1 3r  ^o 

O ifi-  ’ O . Q - O 1J© 

r 3^  r a in  3<f>  T 

5 ,a-  « 

K D.lrf]  ; 


B = B°  + h *L.  B° 

v 3r°  r 


E-.  = r ££  B°  + 21-  ^ B°  +•  £— 

3r°  r r°  3i^  r0sin’^  3<p°  ^ 


0I>  „o 

D. 


B_  = r s 

<P 


lnl^S£_  5°  4-  h sirnf-ag  Bc  + r ain.ifr  3£_  Bo 

“ O _ O J#S-  O p o ® 

r 3£-  r sini?-  9f  Y 


3r 


c r 


54- 


and  for  t ho  vector  potential 

_ 3r°  ko  , o d'J?  4o  o . -o  d®''  4o 
- TT  A_  + r A^.  + r sin^  x*-  A,. 


dr 


Tr~  A<p 


_ 1 dr  "4o  r 3>  4o  . r sintf-  dg>~  , K oo  ^ 

a t tp  K * — ay  t*-1-  — ? — ■&  A?  ( ■ -22) 


= 


dr°  Ao 


O ==,  O O . JJ  C „ 

\ „ 8 A°  . r sinig,  d»  . c 


l«p  r sini>d<p  Ar  ‘ r sin\$"d<p  A r slnS5  d?  A*p 

In  the  Lagrangian  formulation  of  fluid  motion  we  are  ocnoarned 
with  the  trajectories  of  the  fluid  particles. 


r - r(r',  t ) (5.23) 

<VM  /*»** 

Given  these  three  relationships  we  can  compute  the  partial  dev- 
iatives  appearing  in  the  preceding  formulas.  As  a rule,  however, 
we  do  not  start  from  a set  of  trajectories  but  from  a velocity 
field,  v,(r,  t),  corresponding  more  closely  tc  the  Euierian  point 
of  view.  If  the  fluid  motion  is  stationary  we  have  dv/dt  = 0 
and  the  trajectories  (5.23)  are  the  solutions  of  the  differential 
equations 

djy'dt  = £<*>  X0)  (5.24) 

where  It  is  assumed  that  the  Initial  positions  of  the  particles 
may  appear  as  parameters. 

We  may  Inquire  into  cases  in  which  (5.24)  can  be  inte- 
grated by  general  methods.  One  such  case  is 


dx/dt  = v (x,  x°)  , dy/dt  = v (y,  x°  ) 

y 

dz/dt  = v ( z , x?  ) 

S X 

The  integration  can  be  carried  out  at  once  in  the  form 


( 5.25) 


t = / dx/vx  = / dy/v^  - / dz,/vz 

J 


( 5.2S) 
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whence  the  traj oc tcries  (5.23)  follow  by  inversion  of  the 
functional  relationship  between  t ana  x.  etc.  Now  the  dif- 
ferential equations  (5.25)  are  purely  algebraic  relationships 
with  respect  to  the  dependence  of  the  upon  the  x^j  hence  these 
equations  and  their  integrals  may  be  applied  to  curvi 3 ^noar  co- 
ordinates without  further  complications.  Among  the  solutions  are 
in  particular  heiical  trajectories  such  a3  a cylindrical  helix  or 
a helix  winding  along  a cone.  Some  of  these  motions  will  appear 
in  our  later  analysis. 

In  order  to  get  a clearer  though  somewhat  elementary  con- 
ception of  amplif icatory  processes  we  next  discuss  the  effects 
of  induction  for  a field  trial  is  homogeneous  over  some  region  of 
space.  By  (5.17)  the  deformation  of  the  field  depends  only  on 
the  components  dx^/3x^  of  the  strain  tensor.  (This  terminology 
agrees  with  t-be  conventional  definition  of  strain*  the  strains 
are  here  finite  for  finite  times.)  If,  as  usual,  we  disregard  a 
pure  rotation  of  the  fluid  os  if  solidified  (in  'which  case  it  is 
readily  shown  that  the  field  rotates  wi th  the  fluid)  we  are  left 
with  a symmetrical  strain  tensor  which  may  be  decomposed  into  a 
pure  dilatation  (expjmsion  or  compression)  and  a pure  shear. 
Consider  first  a pure  shearing  strain*  We  may  in  this  case 
assume  the  fluid  to  be  incompressible.  If  we  let  the  original 
field  be  in  the  z-directicn  and  the  fluid  motion  in  the  x-direction 
(5.17)  reduces  to 

° (5.27) 


Bx  = Bz  Ox/cK  °) 


» B — 0 « 
* y * 


The  field  energy  is 

1/2  B2  = 1/2  (B°)2[l  +■  ( dx/Bz°  fj 


( 5.28) 


-56- 


The  field  increases  linearly  with  the  strain.  We  see  that  ampli- 
fication of  a field  occurs  whenever  a velocity  shear  field  is 
superposed  perpendicularly  upon  an  existing  magnetic  field,  i'he 
amplification  is  linear  in  time-  for  a stationary  motion  as  appears 
if  w<=*  write  the  first  of  (5.27)  in  the  form 

Bx  = Br  c (9V*Z) 

This  is  illustrated  in  Pig.  1 which  shows  on  the  left  the  velo- 
city profile  and  on  the  right  the  ” stretching"  of  the  original 
(dashed)  magnetic  lines  of  force  in  the  x-direction.  In  principle, 
an  ampiif icat ory  process  of  this  kind  may  be  continued  indefinitely, 
especially  If  the  velocity  field  is  circular  (Sec.  6)  rather  than 
along  straight  lines. 


Figure  1 

It  is  interesting  to  note  that  this  simple  process  by 
(5.28)  always  gives  ar:  increase  of  the  field  energy.  In  order  to 
see  what  happens  to  the  energy  in  more  general  cases,  let  us  return 


( — 


to  (5.17)  and  transform  the  stress  tensor  to  diagonal  form  by  a 
! suitable  rotation  of  the  cartesian  axes  to  which  the  refer 

(we  leave  the  axes  to  which  the  x°  refer  unchanged).  Let  B.,  now 
refer  to  the  components  of  B in  the  new  axes  (the  B?  remaining 
unchanged) . If  £ ^ a~e  the  elements  of  the  diagonalized  stress 
tensor  the  energy  becomes 

1/2  5Z  b5  = J-/2  (5.29) 

1 “ i 1 1 

whereas  the  condition  of  incompressibility  (pure  shear)  merely 
requires  £.  ^ = G.  Now  (5.29)  can  correspond  to  e decrease  or 

an  increase  of  the  magnotic  energy,  depending  on  the  magnitude 
of  the  £ but  for  sufficiently  large  strains  there  will  be  an 
increase.  This  tendency  of  the  magnetic  energy  toward  increase 
is  in  the  long  run  offset  by  the  action  of  the  p order* emotive 
forces  (2.17)  when  they  grow  large  as  the  field  increases;  By  a 
» well-known  principle  these  forces  act  always  in  such  a way  that 

they  In  turn  tend  to  decrease  the  magnetic  field  energy.  (As 
pointed  out  already,  the  dynamo  theory  is  not  based  upon  a trend 
toward  statistical  equilibrium  hereby  implied,,  but  on  certain 
general  dynamical  principles  which  will  be  discussed  later.) 

Consider  the  case  of  a pure  dilatation  (expansion  or  con- 
traction) in  the  absence  of  shear*  Bet  be  in  the  z-direction 
and  assume  for  simplicity  that  the  stress  tensor  1b  diagonal 
along  the  cartesian  axes  given.  Then  (5.17)  reduces  to 

♦ 

Bz  =*  (dz/dz°)  (5.30) 

« 

We  may  exemplify  this  by  a gas  cloud  containing  a homogeneous 
magnetic  field  (for  instance  a 


cloud 


of  ionized  gas  shot  out  by 
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th.a  sun)  expanding  into  a vacuum*  ^ in  o total  magnetic  energy  of 
the  homogeneous  cloud  in  proportional  to  B^/p,  If  we  let  the 
cloud  expand  laterally,  constraining  the  particles  so  that  they 
cannot  move  in  the  z -direction;,  the  magnetic  energy  diminishes 
proportional  to  (jf  % if  ttie  cloud  in  pi-event sd  from  lateral  expansion 
and  expands  only  along  the  z-axis,  the  magnetic  energy  goes  as 

Q -1 

~ , thus  it  would  increase  on  expansion}  if  the  cloud  expands 

laterally  and  contracts  longitudinally  (as  suggested  by  the  dir©ction 
of  the  Maxwellian  stresses)  in  such  a way  that  the  density  remains 
unchanged,  the  magnetic  energy  changes  as  (Az)"  where  Az  is  the 
change  in  extension  in  the  z -dir ection*  For  a gas,  the  internal 
thermal  energy  changes  on  expansion  as  f where  Y'  is  the  ratio 
of  the  specific  heats*  under  suitable  constraints  the  magnetic 
energy  may  thus  incraasa  at  the  expense  of  the  thermal  energy.  We 
might  remark  that  for  a volume  with  a homogeneous  magnetic  field 
the  ponderomotive  forces  (2.17)  vanish  on  the  inside}  they  appear 
as  stresses  on  the  boundary  of  the  volume  and  are  transmitted  to 
the  inside  through  hydrostatic  pres.;  _ces. 
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